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SuMMARY: Values of the pressure coefficient C, between + 1 and - 1 can be 
measured directly in incompressible flow by comparison with the known 
pressures on a circular cylinder placed in the tunnel working section. The 
development of the method and the apparatus required are described. 
The method gives more consistent and accurate readings than can easily be 
obtained with normal manometer measurements, particularly when the 
pressure differences are small and the tunnel speed not very steady. Since 
the coefficients depend on Reynolds number only, less precise measurement 
and control of tunnel speed are required, unless the coefficients to be 
measured are subject to critical scale effects. Possible ways of extending 
the range of coefficients which can be measured are suggested. 


1. Introduction 


In almost all aerodynamic experiments which involve pressure measurements, 
a pressure ratio, or a ratio of pressure differences, is required, rather than the actual 
value of the pressure. In high-speed work the quantity required is sometimes the 
ratio of the pressure on the model to the stagnation pressure, but in low-speed work 
pressure coefficients—ratios of two pressure differences—are almost always used. 


In view of this it seems surprising that methods for measuring these ratios 
directly are so rarely used. The only instruments for this purpose in general use 
seem to be the various types of Machmeter; these are normally used only for setting 
or measuring the free stream Mach numbers of wind tunnels or aeroplanes, not for 
routine pressure measurements. The reason why such pressure ratio instruments are 
so little used clearly lies in the simplicity and convenience of liquid manometers, 
particularly of multiple manometers for rapid reading of pressures at many different 
points. Such manometers have, however, some disadvantages which are worth 
remembering. With an ordinary multi-tube manometer, each result involves three 
meniscus readings, each of which contains an error. This error arises from the 
reading error and from the effects of surface tension, clamping the tubes, before 
reading, and errors in the tube position in the case of an inclined manometer. The 
accuracy depends largely on the magnitude of the pressure differences to be 
measured and becomes poor if these are small. 
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DEGREES FROM STAGNATION POINT 


Fig. 1. 


Pressure distributions round a circular cylinder (Ref. 3). 


When small pressure differences are to be measured, when great accuracy is 
required or when a single Pitot or static tube is traversed, it is necessary to use a 
single manometer, usually of the type in which one limb is raised on a micrometer 
screw. The tunnel speed usually varies slightly with time, so it is often necessary 
to have a second observer measuring the tunnel dynamic pressure, which is to be 
used in reducing the pressure measurements to coefficient form. Even then, rapid 
fluctuations in tunnel speed appear as random errors in the results. The actual 
coefficients may be much more steady, since they would be affected only by 
fluctuations of flow direction, and not by speed fluctuations. 
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SPOILER 


CYLINDER 
\ 12 IN. 
FLOW 
A LAMINAR FLOW PLATE 
WORKING SECTION 


Fig. 2. 


The tunnel, showing the cylinder position for the experiment on the flow 
past spoilers attached to a flat plate. 


Thus it appears that a pressure ratio instrument would have two main 
advantages, namely : — 


(i) Accuracy and ease of operation under unsteady conditions. Only one 
operator is required as against two for accurate work with normal manometers. 
With such manometers high accuracy may be obtained if the tunnel conditions are 
steady enough to allow the use of very sensitive gauges; these gauges, however, 
require skilled operators. A null indicator gives the same sensitivity more simply, 
and should give high accuracy with less experimental difficulty and under less 
perfect conditions. 


(ii) Independence of tunnel speed. The pressure coefficient, at incompressible 
flow speeds, is a function of Reynolds number only, so unless there are critical 
Reynolds number effects only a rough estimate of tunnel speed is needed, and the 
speed need not be kept accurately constant. It should be noted, however, that 
the lack of sensitivity to tunnel speed and hence to speed fluctuations is due to the 
fact that these have the same effect on the model and on the instrument. The same 
will not normally be true of fluctuations in flow direction, and the presence of such 
fluctuations will limit the sensitivity of the instrument. 


In the experiment for which the instrument was developed, the pressure rise 
ahead of a spoiler in the laminar boundary layer of a flat plate was being measured. 
This involved measuring coefficients of the order of 0-1 at speeds down to 15 ft./sec. 
So to get 5 per cent. accuracy the pressure measurement would have to be accurate 
to 3x 10~* in. of water. At first a Chattock gauge was used, but the experimental 
scatter was large, the results were clearly erratic and could not be consistently 
repeated. - It was therefore decided to replace the Chattock by some instrument 
which would be easier to use and also more sensitive. The only way of achieving 
this seemed to be to use a pressure ratio instrument combined with a null indicator. 
The instrument developed and described here enabled accurate measurements to be 
made, and with reasonable ease. The principle involved may have wider applica- 
tions and some possibilities for its development are discussed. 
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NOTATION 
C, a pressure coefficient, =(P—P,)/(H—P,) 
P local static pressure at the point under consideration 
P, free stream static pressure 
H _ free stream total pressure 
6 the angle between the cylinder radius through the pressure hole and 


that through the zero C, position 


na factor depending on the cylinder Reynolds number 
(see Section 4.1(3) ) 


ma factor depending on the thickness of the wake produced by the 
model (see Section 4.1(4)) 


2. The Circular Cylinder as a Pressure Coefficient Meter 


Pressure distributions round circular cylinders have been measured by Schiller 
and Linke"), and by Fage and Falkner®-*’, some of whose results are’ shown in 
Fig. 1. These show that the pressure coefficient varies from +1 at the stagnation 
point to about —1:3 just ahead of the laminar separation. Suppose that such a 
cylinder, with one pressure hole in its surface, is placed in the working section of the 
tunnel in which tests are being made on a model. The speed of flow past the 
cylinder will automatically be proportional to that past the model (subject to 
certain reservations noted later). So if the cylinder is rotated so that the pressure 
at its pressure hole is the same as that at a point on the model, the angle of rotation 
(from the zero C, position) will be a measure of the pressure coefficient. This is 
the basis of the instrument to be described. 


QUADRANT A 


CYLINDER AND 
CONNECTING TUBE 


ZERO SETTING SCREW C 


MICROMETER SCREW B 


Fig. 3. 
The cylinder and adjusting screws. 
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MEASUREMENT OF PRESSURE COEFFICIENTS 


TO CYLINDER 
TO PRESSURE HOLE 


CONNECTING TAP 


U-TUBE OR NULL INDICATOR 


6) 


Connections to the null indicator. 


3. Apparatus and Experimental Details 


3.1. THE CYLINDER 


The apparatus was developed in a 12 in. x 8 in. open-return tunnel which gives 
a maximum speed of 90 ft./sec. The instrument isan approximately circular and 
cylindrical length of copper gas-pipe, } in. diameter, spanning the 8 in. tunnel width 
and rotating in two holes in the wooden tunnel walls. It is placed at the rear of the 
working section near the centre line of the tunnel, as shown in Fig. 2. The cylinder 
has one pressure hole, about 0-01 in. diameter, in the surface. One end of the 
cylinder is fixed to the quadrant A as shown in Fig. 3, and it can be rotated about 
its axis by the micrometer screw B, or by the zero-setting screw C. The micrometer 
screw is 4§ in. from the axis, so a travel of 0:1 in. corresponds to a rotation 
of 1-27°. 


Initially the cylinder was set relative to the quadrant so that, with screw B 
at zero and C at about its mid-position, the pressure at the hole was approximately 
equal to the tunnel static pressure. It was then fixed to the quadrant and further 
adjustments were made entirely by the two screws. 


In use, the zero-setting screw C is adjusted until the pressure at the hole is 
equal to the tunnel static or datum pressure, with screw B at zero. For this purpose 
the hollow cylinder is connected to one side of a null pressure indicator; the other 
side is connected to the datum pressure hole. Next this datum pressure hole is 
disconnected, and the indicator is connected to the pressure to be measured. The 
micrometer screw B is adjusted until the cylinder pressure is equal to this, and 
the micrometer reading then gives a measure of the pressure coefficient required. The 
effect of the choice of the datum pressure is discussed in Section 4.2. 
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Fig. 5. 
The null indicator. Sectional view. Half size. 


3.2. THE NULL INDICATOR 


The type of indicator used for detecting a pressure difference will depend 
mainly on the magnitude of the pressure differences involved. For most ordinary 
low-speed work a simple inclined U-tube, with a tap to connect the limbs, as shown 
in Fig. 4, would probably be suitable. For the experiment considered here, however, 
a much more sensitive device was required, and Fig. 5 shows a sectional view of the 
type used. Essentially this is also a U-tube. The two limbs are brass capsules 
1} in. diameter and the connecting channel, about } in. x 7; in. in cross section, has 
its roof and one side made of Perspex. A beam of light from a projector lamp 
illuminates particles in the fluid in this channel, and these are viewed through a low 
power microscope as shown. The area ratio is such that a pressure difference of 
0:001 in. of water gives a mean particle movement of about 0:06 in., and the 
corresponding travel in the microscope field is about 2 in. It is therefore easy to 
detect a pressure difference of 10~* in. of water. The sensitivity, in this sense, could 
be increased almost indefinitely by increasing the area ratio, or the power of the 
microscope, but in use the sensitivity is limited by the fluctuations in the stream. 


The fluid used is water with the addition of two or three drops of milk to 
provide suitable particles. This needs changing every few days as the particles 
coagulate. Possibly a more stable material could be found, but the milk is the 
best of those tried. The capsules have screw tops so that the fluid can be 
easily replaced. 


One side of the instrument is connected to the cylinder hole and the other to 
the pressure hole with which this is to be compared. A tap, as in Fig. 4, is arranged 
so that the two sides may also be connected directly together. The cylinder is 
rotated until opening or closing the tap causes no appreciable movement of the 
particles. This is a trial and error process which is, however, very rapidly con- 
vergent, particularly when the required cylinder position is roughly known as, for 
example, when taking readings at consecutive holes along the plate. The main 
practical difficulty is due to the fluctuations. It is best not to try to observe the 
mean direction of motion when the tap is closed, but rather to check that there is 
little movement on opening the tap. 
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Obtaining a reading takes about as long as with a Chattock gauge, owing to the 
trial and error process and the fluctuations. The apparatus, however, follows 
pressure changes much more rapidly than a normal Chattock, as is evident from 
the fact that the null indicator is affected by rapid fluctuations. This is probably 
due mainly to the fact that the connecting channel of the null indicator can be 
much smaller in area than the Chattock meniscus, so that less air has to flow through 
the pressure hole to cause a given movement. 


4. The Calibration of the Cylinder, and the Assumptions Involved 


4.1. THE CALIBRATION MEASUREMENTS REQUIRED 
There are four points to check in calibrating the cylinder, namely :— 


(1) The position of zero pressure coefficient. This depends on tunnel speed 
(i.e. on cylinder Reynolds number) and on the model in the tunnel, since this may 
alter the flow direction at the cylinder. To allow for this the cylinder is adjusted by 
the zero-setting screw until the pressure is equal to the datum static pressure. This 
is normally done at each run of the tunnel, although this may not be necessary if 
the speed and model are the same in the next run. 


(2) The variation of pressure coefficient with micrometer reading. This will 
probably be linear if the required range is small; otherwise a complete calibration 
is necessary. 


(3) The variation of micrometer reading with Reynolds number, at a given 
pressure coefficient. 


The pressure distribution round the front part of a circular cylinder varies very 
little with Reynolds number, so long as this is always below, or always above, that 
of the change from laminar to turbulent separation. The cylinder size used should 
be such as to avoid the critical range of Reynolds number. Normally, Reynolds - 
numbers below the critical would be used, except perhaps in a compressed air 
tunnel. For example in a normal tunnel at 300 ft./sec. and with a cylinder diameter 
of 1 in., the Reynolds number would be about 160,000, which would probably still 
be below the critical value. With the cylinder diameter of 3 in. and the tunnel 
speeds used in the present experiment, the Reynolds number ranges from 3,000 
to 18,000. 


Points (2) and (3) are properties of a given cylinder and tunnel only, and to a 
first order are independent of the model under test. So it would not be too difficult 
to obtain a complete set of readings giving C, as a function of micrometer reading 
and Reynolds number. This is not necessary, however. The variation with 
Reynolds number is such that it can be assumed, for C, between +1 and —1 at 
least, that the pressure distribution has the same shape at all Reynolds numbers, the 
sensitivity (rate of variation of C, with 9) alone varying. That is, it can be assumed 
that C, is a function of n6 or n x the micrometer reading, where n is a factor depend- 
ing on the Reynolds number. It is then only necessary to measure n, which 
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involves one measurement only at each speed. This can be the reading at the 
highest C,, positive or negative, in the required range, which can be done accurately 
even at low speeds. The shape of the distribution, (2), can be determined with good 
accuracy at the highest speed. 


Alternatively, the factor n could be determined by measuring the angle between 
the two positions of zero C,, one on each side of the front stagnation point. No 
manometer readings would then be needed except at the highest tunnel speed, but 
a device for measuring the angle of rotation over the range of about 70° would 
be required. 


Thus the calibration measurements could all be done at the highest tunnel speed, 
and so the accuracy of the method is not limited to that obtainable with normal 
manometers, even though these have to be used for the calibration. 


(4) Variation of sensitivity with the configuration of the model. There are 
two possible effects namely: — 


~ 


(a) If the model produces a wake, the sensitivity will increase, since the 
velocity at the cylinder will be raised. 


Once again the shape of the curve of C, against 6 should be unchanged, 
so C, will be a function of m#, where m depends on the thickness of the 
wake. This will depend on the model and possibly also on its Reynolds 
number, as for example with a wing stalled at low speeds but not at high 
speeds. Again, if m has to be measured at low speeds this could be done 
most accurately by measuring the angle between the two positions of 
zero C,. 


This dependence on the wake produced by the model could be 
eliminated by placing the cylinder ahead of the model, at the front instead 
of the rear of the working section. The front position would thus be 
better, so long as no trouble arises from the turbulent wake of the cylinder. 


(b) If the cylinder is too near the model, the model will produce a velocity 
gradient affecting the cylinder pressure distribution. The sensitivity will 
be affected differently on the two sides of the stagnation point, so the 
change cannot be measured by measuring the angle between the zero C, 
positions. In practice the cylinder should be made small and placed 
sufficiently far from the model to ensure that the effect will be negligible. 


These effects will probably be most serious if the model is an aerofoil at 
incidence. The flow direction at the cylinder will vary with the aerofoil incidence, 
so the zero C, position must be readjusted after any change of incidence. The most 
serious effect, however, is that the circulation round the aerofoil produces a velocity 
gradient at the cylinder. In the present tunnel, if the 3 in. diameter cylinder 
were placed 12 in. from the centre of an aerofoil of 4 in. chord, at a lift coefficient of 
unity, the variation of velocity from one side of the cylinder to the other would be 
roughly 0:001 of the stream velocity. Thus the effect can be made fairly small. 
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4.2. THE METHOD OF CALIBRATION 

The cylinder is first adjusted so that the pressure at its pressure hole is equal to 
the datum static pressure, by means of the zero-setting screw. The tunnel dynamic 
pressure is then measured on a Chattock or other manometer. Then the Chattock is 
set to a definite fraction, for example {, of the reading obtained, and the cylinder is 
connected to it in place of the Pitot tube. The micrometer screw is then adjusted to 
balance the Chattock, and the reading obtained is that corresponding to a pressure 
coefficient of 0-2. Finally the Pitot tube is replaced to check the dynamic pressure. 


It should be noted that the pressure coefficient measured is (P—P,)/(H —P.), 
where P, is the datum pressure used for setting the cylinder when taking measure- 
ments on the model, and P, is the datum used in the calibration. The method 
should give an accurate measure of this quantity. It is an accurate measure of C, 
only if P, and P, are both equal to the true free stream static pressure, and this will 
depend on the design and location of the static tubes or holes in the usual way. 


5. Results 


5.1. CALIBRATION 

The range of pressure coefficients required for the experiment on the flow 
ahead of spoilers was small, and so calibration was carried out over the range of 
C, from 0 to +0-2 only. The results at various speeds given in Table I show that 
the variation with micrometer reading was linear and that the sensitivity was 
independent of Reynolds number in this range (3,000 to 12,000). The factor n is 
thus constant and may be taken as unity. 


TABLE I 
READINGS FOR C,=0'1 AND 0:2 AT VARIOUS SPEEDS, WITH A } IN. SPOILER IN POSITION 
(in. water) Micrometer Reading 
Reading Reading C,=02 C,=01 C,=02 
0:0685 0:0682 0-167 0:340 0°599 0°588 
0°1028 0°1028 0-170 0°330 0°588 0°606 
0°171 0°172 0-167 0°330 0-599 0°606 
0:360 0°361 0°167 0°330 0:599 0-606 
0°627 0°626 . 0:330 | 0-606 0°606 
0-912 0-915 0-167 0:335 | 0-599 0°596 
| | Mean 0-600 
TABLE II 
READINGS AT TOP SPEED AND C,=0:2, WITH VARIOUS SPOILERS 
Spoiler Dynamic Pressure | Cylinder | Micrometer | C, 
Height (in. water) Pressure Reading | Micrometer Reading 
Before After | 
Reading Reading | 
No Spoiler 0°884 0-880 | 0179 0354 | 0-571 
4 in. 0858 | 0852 | 0172 0-348 | 0:576 
fin. 0852 | 0852 | O17] 0-338 0°593 
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fe) 
16 FT/SEC. 
= 4 
0-01 (H-P,) =0°65 x 10° IN. WATER 
WwW 
8 ol NX 
> 24FT,/SEC. 
| 0 IN. WATER \ 
a 


0-2 


0-01 (H-P,) =3°25 x10 IN. WATER \ 
0-2 
+ 
0-3 
io 9 8 7 6 a 4 3 2 1 re) 


INCHES UPSTREAM OF THE SPOILER 
(Leading edge of plate 12 in. upstream) 
“Fig. 6. 
Pressure distributions upstream of a } in. spoiler on a flat plate with laminar boundary layer. 


A set of readings with different spoilers in position was taken, at the highest 
speed and a C, of 0:2. These show that a } in. high spoiler in the 12 in. high tunnel 
increases the sensitivity by about 4 per cent., corresponding to an increase in dis- 
placement thickness of the wake of about 4 in. The factor m thus varies from unity 
with no spoiler to 0:96 with the 4 in. spoiler. The slight discrepancy between the 
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0-1 
55 FT/SEC. 
0:01 (H-R,) =7-8.x 10° INWATER 
4 
0-2 
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8 fe) \ 
80 FT. /SEC. 
0-01 (H-P,)= 16-3 x 10° IN.WATER 
0-2 
0-3 
0-4 
10 9 8 7 6 5 4 3 2 1 te) 


INCHES UPSTREAM OF THE SPOILER 
(Leading edge of plate 12 in. upstream) 


Fig. 7. 
Pressure distributions upstream of a } in. spoiler on a flat plate with laminar boundary layer. 


two readings given in Tables I and II for the 4 in. spoiler is probably due to a 
difference in the incidence of the plate on the two occasions. 


5.2. PRESSURE DISTRIBUTIONS 


Figures 6 and 7 give pressure distributions at various speeds from 16 to 
80 ft./sec., which show that the readings are consistent and the random errors small, 
even at low speeds and small pressure coefficients. The pressure difference corre- 
sponding to a coefficient of 0-01 is shown on each graph, and ranges from 
0-65 x 10-* in. to 16-3 x 10-* in. of water. 
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MODEL VENTURI TUBE 


Fig. 8. 
Methods of increasing the measurable range of C,. 


5.3. GENERAL BEHAVIOUR OF THE APPARATUS 


It is perhaps worth noting that the cylinder part of the apparatus was assembled 
in a few hours, that no modifications or development were needed, and that no 
trouble has been experienced in routine experimental use. The development of the 
null indicator proved more difficult, but since the design shown in Fig. 5 was made, 
and a suitable fluid found, no trouble has occurred. Despite initial difficulties with 
the null indicator, the complete equipment was in routine use within about a fort- 
night of the initial decision to try the idea. 


6. Further Possibilities for Development 


The arrangement described could be used over a range of C, from +1 to 
about —1. For still larger negative coefficients there seem to be two possible ways 
of extending this range, namely : — 


(i) An aerofoil, with a pressure hole at or near the leading edge on the upper 
surface, could be used instead of the cylinder. This should give large 
negative coefficients, but more care would be needed in calibrating to find 
any Reynolds number effects. In particular, the simplifying assumptions 
made in Section 4.1 in calibrating the cylinder at varying Reynolds num- 
bers would probably not be applicable. 


(ii) A very small auxiliary tunnel could be used, driven by the main tunnel 
fan but with a more efficient diffuser giving a higher speed. This could 
be separate from the main tunnel, except for the end of the diffuser, or 
alternatively it could be in the form of a venturi tube held in the tunnel, as 
Fig. 8 shows. For example, if the small tunnel ran at twice the speed of 
the main tunnel, a circular cylinder placed in it would give a range of C,, 
based on main tunnel speed, from +1 to —7, which should be adequate for 
most uses. Further advantages would be that the sensitivity should depend 
less on the wake produced by the model in the tunnel and that the effect 
of fluctuations in flow direction may be reduced. 


The method could probably be made automatically self-balancing with little 
difficulty, by using some type of null indicator to control a motor rotating the 
cylinder. 
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7. Conclusions 


Pressure coefficients ranging from +1 to —1 can be measured by comparing 
the pressure on the model with that on a circular cylinder placed at the front or rear 
of the working section. The method has two main advantages : —— 


(i) It gives more consistent readings than normal manometer measurements, 


particularly for small pressure differences or in unsteady conditions. The 
instrument is insensitive to fluctuations in tunnel speed, although it is 
affected by fluctuations in flow direction. The calibration needed may 
often be done entirely at high speeds, so the accuracy is practically inde- 
pendent of manometer measurements. Being a null method it can be used 
over a wide range of pressure differences while a very high sensitivity 
is retained. 


d 
0 (ii) Since it measures coefficients directly, the reading usually varies only 
e ; slowly with speed, so no precise measurement or regulation of tunnel 
speed is needed. 
h The possible range of coefficients might be greatly increased by using an 
- ; aerofoil instead of the cylinder, or by placing the cylinder in a small auxiliary tunnel 
j running at higher speed. 
} The apparatus is simple and convenient, requires little development and has 
given very little trouble in use. 
O 
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Compressible Subsonic Flow in 
Two-dimensional Channels 


Part II: The Application of the Theory to Problems of 
Channel Flow* 


L. C. WOODS, D.Phil., A.F.R.Ae.S. 


(Department of Applied Mathematics, University of Sydney) 


SuMMARY: The theory developed in Part I* is applied to a number of 
problems of aeronautical interest, the most important of which is to the 
setting of “streamline” walls about a symmetrical aerofoil placed in the 
centre of a channel. It is shown how the position of the streamline wall 
can be deduced from the (experimentally determined) position of a constant 
pressure wall. This theory is applicable to symmetrical aerofoils of any 
given shape, and makes allowance for the presence of the aerofoil’s wake. 
To illustrate the theory, and to test it by an extreme example, the flow is 
calculated about a circular cylinder, with a diameter about half the tunnel 
height, for both straight and constant pressure walls. The solid blockage is 
calculated in each case and compared with the standard first order theory. 
For this extreme example the standard theory fails badly for straight walls, 
but is reasonably accurate for constant pressure walls. 


NotaTION (Additional to that defined in Part I*) 
gd x = 2tan~' (tanh 4x), the Gudermannian function 
R _ radius of curvature of a solid boundary 
A __ cross-sectional area of an aerofoil 
blockage factor, (q—4q:)/q: 
o = y;/y,, On the walls 
c aerofoil chord 


Cy drag coefficient 


” 


Suffixes “p” and “f” denote values appropriate to “constant pressure walls 
and “streamline walls,” respectively. 


7. Summary of Equations 
Let (g, @) be the velocity vector in polar co-ordinates, 9, Y be the velocity 
potential and stream functions, lim (g, #)=(U, 0) (upstream), lim (g, #)=(V, 7), 
o> -& 


*Part I of this paper was published in August 1955 (pp. 205-220). 
Received June 1955. 
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FLOW IN CHANNELS 


M,, be the Mach number at ¢=-0, Pol Po, Where 
Poo» Po are the fluid densities at ¢= — 0O and a stagnation point respectively, 


e=sinh—'(6,,/M,.), ‘ ‘ (1) 


and is defined by = sinh cosech (© + ¢). ‘ (2) 


q 
U 


Then it has been shown in Part I that the flow through a given channel can be 
calculated from 


co 


1 1 1 U 
(3) 
U_ 
g=-& 


where f=0+i6, w=o+im,y, H is the channel width at ¢= —00, and are 
the values of 6 on the walls Y=0, Y=h/m,, respectively. 


It has also been shown that the physical or z(=x+iy)-plane is related to the 


w-plane by 


_ 
Udz= 2B. e’ dw dw, (6) 


and that the channel width K downstream at infinity is given by 


V 


If the wall Y=h/m,, is a constant pressure wall, it has been shown in 
Part I that 


1 
f= 6, cosech ah (o—w)do. . (8) 
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Also, for channel design the equations 


== (g—w)—§ } 


have been established. 
One final result required in this paper is that when @,(@) is differentiable 


we can write 


where R(=ds/d0) is the radius of curvature of the wall, and q, is the velocity at 
the wall: a similar result holds for d6,,. 


8. Boundary Conditions Prescribed as Simple Step-Functions 


There are many problems in which the boundary values are specified as simple 
step-functions. For example, in the direct problem of calculating the flow through 
a given channel the walls may be composed of a number of straight sections; in 
this case the direction of flow, 6, is a step-function on the walls; Fig. 1 shows such 
a case. A similar case is shown in Fig. 2, except that here one wall is at constant 
pressure. 


z y FREE SURFACE 
Vv a \= 


¢ 
Fig. 1. Fig. 2. 
Channel with rectilinear walls. Flow of a stream up a step. 


Similarly, in the indirect or design problem it may be considered desirable to 
prescribe constant pressure over a number of intervals, for this avoids adverse 
pressure gradients except those occurring at points of pressure discontinuity, and 
flow separation at these points could be prevented by some form of boundary layer 
control. Lighthill”) has used this principle to design wind tunnel bends; Fig. 3 
shows such a bend, obtained by prescribing q=U on E_..CE,, E-.~A and 
BE,,, and q=V on AB. 
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8.1. FLOW THROUGH A CHANNEL WITH RECTILINEAR WALLS 


The following theory is not strictly applicable to compressible flow, as with 
flow past a sharp projecting corner the local Mach number exceeds unity in the 
neighbourhood of the corner; however, outside such neighbourhoods the theory 
will yield reasonably accurate results. 


— that 6, is constant except for discontinuities of amount 7; at ¢=4,, 


j=1, 2,...m, and that 4, is likewise constant except for jumps 7; at ¢=9,, 
i=1, 2,...n, then the Riemann-Stieltjes integrals in (3) and (4) degenerate to give 
f= log 7; log sinh w)+ + log cosh — (9;—w), 
and h log = 


Hence from (2), (5) and (6) the required solution is contained in 


Uz= 2B 


Poe | dw 


| dw— | aces 


=I exp (—7#;/HUB.) 


q=U 


and A(w)= ()’ cosh (¢;- [ sinh (6,-w) 


where the ratio U/V is determined from (7). The constants ¢;, ¢; are calculated 
from (12) and the known corresponding values of z. Notice the simplification that 
occurs when £,,=1, that is when the flow is incompressible. 


4; 


The usual method of solving these problems is by means of the Schwarz- 
Christoffel mapping theorem, which yields the same results, but with more labour. 


8.2. DESIGN OF A CHANNEL WITH WALL PRESSURES SECTIONALLY CONSTANT 


The calculation is similar to that given in Section 8.1. From (9) and (10), 


(13) 


2HUB... 


1 


where ©; are the jumps in 2), at ¢=49,, and {2; are the jumps in ©, at o=4,. 
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As a simple example, suppose 2, jumps from 0 to o at ¢= —a, and from o to 0 
at ¢=a, and that 2,,=0 for all 9; then, from (13) and (14), 


sinh 


As indicated at the beginning of Section 8, these equations can be used to design an 
efficient bend in a channel; « and a are chosen to give the desired value of + 
(equation (16) )}—the angle of the bend—then (15) yields the required wall shape. 


Fig. 3. 
Design of a right-angled bend. 


For example on ¥=0, |¢|<a, equations (15), (1) and (2) yield the intrinsic 
equation 


sinh {t (Vs+a)/2HUB,} 
| sinh (= (Vs—a)/2HUB,.} 


where o=sinh- ay) —sinh-* 


and s is distance measured along the wall. 


The incompressible flow form of these equations was obtained by Lighthill"’ 
by means of the Schwarz-Christoffel theorem. Fig. 3 gives some idea of the type 
of bend achieved. 


8.3. MIXED BOUNDARY CONDITIONS: FLOW OF A STREAM UP A STEP 


Equation (8) can be reduced similarly to give results corresponding to those of 
Sections 8.1 and 8.2; we omit this calculation but consider the simple example 
shown in Fig. 3, the flow of a stream over a right-angled step. 
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FLOW IN CHANNELS 


From equation (8) it is readily found that the intrinsic equation of the free 
surface is 


(Us +a) n(Us—a) 


1 


where a is defined in the figure, and gd x (=2 tan~' (tanh 4x) ) is the Gudermannian 
function. 


The co-ordinates of the free surface are 


z= | cos 6,ds, y= sin 6,ds, (18) 
0 0 


and eliminating s from these, by (17), gives 


(sin? 0’ + k? cos? &)! + k (sin? 0, + k? cos? 6,)! —k 
(sin? 6’ + k? cos* &)!— k (sin? 0, + k? cos? +k 


x(0,)= tog 


2HkB.. 


where =gd [(xa/(2HU§8,,)], k=sinh [a/(2HUB,,)], and F is the elliptic integral 
of the first type. If A is the total depth of the step, then from (20) 


4HkB,. 


which gives the relation between A and a. 


This example shows how soon difficult algebra appears, even in the simplest 
cases. For curved boundaries, even if simple functional relations 6,=9, (¢) were 
assumed, the algebra would be excessive; there is no real alternative to the numerical 
methods given in the next section. 


9. Curved Boundaries 


It has been shown in Section 4.1 of Part I that with given curved walls two 
integral equations must be solved in order to determine the flow; if one wall is 
straight or at constant pressure a single integral equation suffices, and for simplicity 
of exposition it is this case which is considered here. 
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J(y)= | log tanh xdx 


9 


y J (y)—y log y y 
0 0 0 0°50 —0°86000 +426 
0:05 —0-05001 0°55 —0°89656 +375 
0-10 —0°10011 0°60 —0°92937 + 332 
0-15 —0°15037 25 0°65 —0-95886 +295 
0:20 —0-20088 — 33 0°70 — 0°98540 +263 
0°25 —0°25172 — 40 | 0°75 1:00931 +234 
0:30 —0°30296 48 | 0-80 ~ 1-03088 +212 
0°35 — 0°35468 36 0°85 — 1:05033 +189 
0:40 — 0°40696 — 61 | 0:90 — 106789 +169 
0°45 — 0°45985 69 0-95 — 1:08376 +154 
0°50 —0°51343 — 74 1:00 — 1:09809 +138 
0°55 —0°56775 — 80 1:05 —1°11104 +124 
0°60 —0°62287 — 86 1:10 — 1°12275 +113 
0°65 —0°67885 — 90 1:15 — 1°13333 +101 
0°70 — 0°73573 — 94 1:20 — 1:14290 
0°75 — 0°79355 — 1:25 —1°15155 + 82 
0°80 —0°85236 — 102 1:30 — 1:15938 + 75 
0°85 —0°91219 — 105 1:35 — 1:16646 + 67 
0:90 —0:97307 — 108 1:40 — 1°17287 + 62 
0°95 — 1:03503 —110 1-45 — 117866 + 55 
1:00 — 109809 1:50 — 1:18390 + 50 
135 — 1:18864 + 45 
1:60 —1-19293 + 41 
1°65 — 1:19681 + 37 
1:70 — 1:20032 + 33 
— 1:20350 + 31 
1:80 — 1:20637 + 26 
1:85 — 1:20898 + 26 
1-90 — 1°21133 + 22 
1:95 — 121346 + 21 
2:00 — 1:21538 + 17 
2°05 121713 + 18 
2°10 — 1:21870 + 14 
245 — 1:22013 + 14 
2°20 — 1:22142 + 12 
2°25 — 1°22259 + 11 
2°30 — 122365 + 10 
2°35 — 1:22461 + 10 
2°40 — 1:22547 + & 
2°45 — 122625 + 7 
2°50 — 1:22696 + 7 
2°59 — 1:22760 + 6 
Note: For y > 3, 2°60 — 1:22818 + 5 
2°65 — 1:22871 + 6 
J (y)=e-29- 9 decimal places). 
2°80 — 1-23000 + 4 
2°85 — 1:23035 + 3 
2:90 — 1:23067 + 3 
2°95 — 1:23096 
3-00 — 1-23122 + 2 
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9.1. THE NUMERICAL SOLUTION OF THE INTEGRAL EQUATION 
(a) One Wall a Constant Pressure Wall 


First we show how to determine the velocity on the given curved wall ¥=0. 
On ¥=0 equation (8) yields 


2, = 


on integrating by parts and denoting the variable of integration by 9’. Now where 
6, is continuous, (11) is applicable; elsewhere suppose @, has simple discon- 
tinuities of amount 7), at ¢= 9, K=1, 2, . . . m; then (21) gives 


oo 


Equation (2), which gives the relation between {2, and g,, shows just how 
complicated this integral equation is. 


In most practical examples 1/R=0 outside some finite range —-2a<¢< 2a, 
say. If this range is divided into n parts thus (2a, $,,_,5%_iy+++-Pji_ys++-P 145 — 2), 
such that 1/(Rq,) can be taken as constant, say 1/(Rq), in the i interval 
(%,,4+ %_,)» With negligible error (excluding the points where @, is discontinuous), 
then we can write 


2,@)=3 | -2 3 log 


tanh — * | 


where J (y)= f log tanh xdx, and [X]; denotes the first difference, or jump in value 
0 


of the bracketed function over. the i” interval. J(y) is tabulated in Table I. 
If Ay=4[4 {= } A lo tanh — 
and B= (x) B.= 
then the formula for 2, can be written concisely 
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TABLE II 
I(Qy)= gd xdx 
v 
y gd xdx y gd xdx 

0 0 
0 0 1000 | 2-5 2:25911 163 
0-1 000500 993 2-6 2-40060 147 
0-2 0-01993 981 2-7 2:54356 135 
03 0-04467 955 | 2:8 2-68787 121 
0-4 0:07896 925 2-9 2:83339 109 
0°5 0-12250 886 | 3-0 2:98000 101 
0-6 0-17490 844 | 31 3-12762 89 
0:7 0:23574 795 | 3-2 3-27613 81 
0:8 030453 749 3:3 3°42545 714 
0-9 0-38081 697 3-57551 67 
1-0 0:46406 649 3-5 3°72624 61 
11 0-55380 599 | 3-6 3°87758 54 
1-2 064953 553 | 3-7 4:02946 49 
1:3 0-75079 507 3:8 4:18183 46 
1:4 0°85712 465 | 3-9 433466 39 
1-5 0:96810 425 4-0 448788 38 
1-6 1:08333 389 41 4:64148 32 
1:7 1:20245 354 42 4:79540 31 
1:8 1-32511 321 43 4:94963 27 
1-9 1-45098 293 44 5-10413 24 
2-0 1-57978 267 4°5 5-25887 23 
21 171125 241 46 5-41384 19 
2:2 184513 218 47 5-56900 19 
2:3 198119 200 4:8 5-72435 16 
2-4 2:11925 180 4:9 5-87986 16 
5-0 603553 12 

ote: or y> (y)= + Ga-ipP 
ie. I (y)=2e-¥+ = y—1-8319312 (to 7 decimal places). 


Besides (22) there is another equation necessary in the numerical calculations; 


oo 


it is the self-evident result | d6, (¢’)=7, where 7 is the total angle of deflection of 


the channel. By (11) and the approximations already introduced this can be written 
1 n 


i=1 k=1 


Equation (23) must be satisfied exactly if the overall deflection angle = is to be 
achieved. This can be done most simply by modifying the values of B, (see the 
example in Section 9.2). 
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Similarly, from (8) and (11), it is found that the slope of the free surface can 
be calculated from 


and I(y)= [ gdxdx. I(y) is given in Table II. 


Since the values of B; in (22) are unknown until 0, has been calculated, it is 
necessary to solve the equation by iteration. If reasonable values of g,/U are 
assumed at ¢=4¢,;, j=1, 2,...n, then 


sU= | do, ‘ ‘ - & 


enables the distance s to be determined as a function of ¢. R(¢) now follows from 
this relationship and the given wall shape. Then B,=[HU8,,/(Rq)], can be obtained 
and used in (22) to yield values of 2, (¢), and hence g,/U from (2). This calcula- 
tion is now repeated until an iteration does not alter q,/U. 


When this has been done the final values of B; can be used directly in (24) 
to give 9,; the co-ordinates of the constant pressure wall then follow from 


Two further points in connection with the application of equations (25) and 
(26) need be mentioned. 


First, near a stagnation point, at ¢=b, say, produced by a jump in @ of r, U/q, 
becomes infinite, when instead of (26) the approximate equation?’ 


must be used to calculate the distance s from the stagnation point to ¢=¢’, where 
the velocity is q’. 
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Secondly, since (%,—7) will in general vanish only at ~=00, it is useful 
to have a formula from which to calculate y at some large value of #, say 9%). 


In this case the second equation of (26) can be replaced by 


Uy=Uy(9,)+ | sin 6,d9. 


Consider, for simplicity, the case =0; then 


Uy (%,)=Uy (00) + sin 9,do 


oo 


% 
=Uy (00) + | 


oo 


provided », is “reasonably” large (see later). Therefore inserting the value of @, 
given by the imaginary part of (8) gives 


y(9,)—y(O)= { sech 5; (9’ —9) dq’ a 


2a 


~ 2a 


Integration by parts gives 


where / is the function already defined, and (5) has been used. Now for x> 4, 
I (x)=2e-*+ 42x — 1-8319 is accurate to four decimal places, and so in the calcula- 
tions it is sufficient to take 


Hi p } 


provided that = (9, +2a)/(2HUB,,.) > 4. 
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Thus in the notation of the first part of this section we can write 


(28) 


(b) One Wall Straight 


Suppose #,=0, then the velocity on the other wall (/=0) can be found from 


m 


where B;, B,, have already been defined, 


and F (= | log sinh xdx. This function is given in Table III. Equation (29) 


follows from (3) and (11) by making the same approximations as in (a). 


In the present case, besides the auxiliary equation (23), it is also necessary 
to satisfy 


43 U_ 
which comes from (4). The ratio U/V can be determined from (7). 


9.2. EXAMPLE: FLOW ABOUT A CIRCULAR CYLINDER IN A CHANNEL 


Consider the incompressible flow about a circular cylinder placed centrally in 
a stream of total width 2H at infinity. As the flow is symmetrical about the 
stagnation streamline, say Y=0, only that half of the flow pattern lying between 
the streamlines Y=0 and y=AU need be considered; this can be regarded as a 
“channel” flow. Compressible flows can be determined by making only minor 
modifications to the following calculations. 


The flow about a circular cylinder may seem to be an unpractical example, but 
it does serve to test the theory severely, for the larger the thickness ratio of the 
body the more slowly does the iterative process described in Section 9.1 converge. 
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TABLE 


F (y)= log sinh x dx | 


0 


y F (y)—y log y | y F (y) i? 
0 0 0 0:50 —0:83966 541 | 
0:05 —0-04999 4 0:55 —0°86962 500 
0-10 —0-09994 8 0:60 —0°89458 466 
0-15 —0-14981 12 0-65 —0-91488 438 | 
0:20 —0-19956 18 0:70 —0-93080 414 
0:25 —0-24913 20 0-75 —0-94258 394 
0:30 —0-29850 25 0-80 ~0-95042 376 
0:35 —0°34762 28 0-85 —0:95450 362 
0-40 — 039646 34 0-90 095496 350 
0:45 —0:44496 37 0:95 —0-95192 337 
0:50 —0-49309 41 | 1-00 —0-94551 329 
0:55 —0:54081 44 1:05 —0:93581 320 
0-60 —0-58809 50 1:10 —0-92291 312 
0°65 —0-63487 52 1-15 — 090689 306 
0:70 —0-68113 57 1:20 —0°88781 300 
0:75 —0°72682 60 1:25 ~0°86573 294 
0:80 —0°77191 64 1:30 —0-84071 291 
0:85 —0°81636 68 1:35 —0-81278 286 
0-90 —0-86013 71 1-40 —0:78199 282 
0-95 —0-90319 14 1-45 — 074838 279 
1-00 —0-94551 79 1:50 —0°71198 277 
1°55 273 
1:60 —0-63091 271 
1-65 —0-58630 270 
1:70 —0-53899 267 | 
1:75 —0-48901 265 
1:80 —0-43638 265 
1:85 —0°38110 262 
1:90 —0-32320 262 i 
1:95 —0:26268 260 
2:00 —0-19956 259 
2-05 ~0-13385 259 
2-10 --0:06555 257 
2-15 +0-00532 258 
2:20 007877 255 | 
2:25 015477 256 
2:30 0:23333 256 
2-35 0-31445 253 
2-40 0:39810 255 
2:45 0-48430 254 
2:50 0-57304 253 
2:55 0°66431 253 j 
2:60 0°75811 253 
2:65 0°85444 253 
Note: For y >3, 2:70 0-95330 251 \ 
2 
F (y)= —y log, 2 + 4y? + 4e- 24 — 
2:85 1:26499 251 
(to 5 decimal places). 5-90 137392 252 
2:95 1:48537 251 
3-00 1-59933 252 
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(a) One Wall a Constant Pressure Wall 


This case is obtained if the cylinder (radius R) is placed centrally in a jet of 
fluid, as shown in Fig. 4. The deflection angle 7 is zero for this case. It is 
assumed that 2a= AU, which yields a ratio of approximately 1/2 for R/H, as in a 
stream of infinite width it is a standard result that R=a/U. 


The range (2a, —2a) is divided into sixteen equal parts, while 9; are taken to 
be at the mid-points of these intervals (see the first column of Table IV). Aj; are 
then simply the first differences of 


= 4 5 
» BAG 


A good initial assumption for a U(¢) is the infinite stream distribution q,/U, 


namely 


From the definition of B; and (5) 


B= (Re), = Ro) 


With the assumption (31), RU =a, and hence 


For the present example @ has two discontinuities (r,), one at each stagnation point, 


and hence (23) yields 
(a0) Al ~ 2B. 


32 


Thus 
as B,=—r,. As the number of intervals into which (2a, —2a) is divided tends to 
infinity, 7, tends to 4x; for the division adopted here 7, is found to be 82°31°. An 
application of (22), in which 2,=log(U/q,) as the flow is incompressible, com- 
pletes the first iteration. Notice that in this example it is unnecessary to use (25) 
during the iterative process, due to the constancy of R. 
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Fig. 4. 
Flow about a circular cylinder. 


The final values of g,/U, obtained after five iterations, are given in column 1 
of Table IV. Column 2 gives the values of sU/(4a), obtained from column 1 and 
equation (25). From the first entry in this column it follows that 


RU /(2a)=4 x 0-4145/x=0-5278, and so R/H=0°5278 as HU =2a. 


Column 3 gives the value of 2, defined in Fig. 4, while column 4 gives the free 
stream value g;=2U sin a. 


Column 5 sets out the solid blockage «, defined by « =(q—4q:)/q:, which agrees 
closely with the first order term used in the standard theory. In the general case 
this first order term can be calculated as follows. 


For slender cylinders 6, ~ dy/dx, ¢~Ux, so on ¥=0 equations (5) and (8) 
yield 


1 
2HB. (x x)dy(x ), 


where c is the “chord” length of the cylinder. Allowing H to tend to infinity gives 


c/2 
1 dy(x’) 


2’=—c/2 
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Now to first order in 2, from (1), (2) and the definition of < it is found that 


e/2 


=—c/2 


e/2 


= | (x’ — x)dy(x’), 


2’=-c/2 
provided that c/H is small. An integration by parts now gives 


TA 
~ 72 (3) 


t= 


where A is the cross-sectional area of the cylinder. 


This theory would appear to be applicable only to conventional aerofoil shapes 
in jets of width large compared to the aerofoil chord. For the circular cylinder 


(8) 


The average value in column 4 of Table IV is 0-0525, and so, even for the extreme 
case of the flow of a relatively narrow jet about a circular cylinder, the first order 
term is still a fair approximation. 


Column 6 of Table IV is obtained directly from equation (24), while columns 7 
and 8 follow from equations (26). The value of y,/H at ¢/(4a)=97/32, was cal- 
culated from equation (28). As a check equation (28) was also used to calculate 
y,/H at »/(4a)=53/32 and ¢/(4a)=33/32, with the satisfactory results 0-0038 and 
0-0399 respectively (cf. the figures marked with an asterisk in column 7). 


In an infinite stream 


. R?*(x-iy) 
¢+iv=U {x+iy+ 
2 
and so on wW=AU, 


which enables column 9 to be calculated. Column 10 is the ratio s=y,/y,, and this 
ratio becomes unity at x/H=1-532, which is in close agreement with the value 
deduced “or an aerofoil of the same chord in Section 10. 


270 The Aeronautical Quarterly 


| 

| 
| 


hat 


33) 


Des 
Jer 


me 
Jer 


7 
al- 
ate 
nd 


his 
ue 


erly 


| 


FLOW IN CHANNELS 


TABLE V 
FLOW ABOUT A CYLINDER WITHIN STRAIGHT WALLS: R/H=0-°4749 

1 2 a 

: Solid 
q sU | \ (infinite 
32 (2 = ) stream) 

0 2:5609 0°3141 90-00° 2:0000 0-2805 
1 2°5551 0°3019 86°50° 1:9963 0°2799 
3 2°5087 0:2772 1-9662 0°2759 
5 2°4142 0-2519 72°18° | 1:9040 0:2680 
7 2°2696 0°2252 64:°53° } 1:8055 0:2570 
9 2:0687 0°1964 56°29° 1-6638 0°2434 
11 1-7988 0:1642 47:04° 1-4637 0°2289 
13 1:4273 0°1255 35-95° 1:1743 0°2154 
15 0°8315 0-0691 19°81° 0:6779 0°2266 

16 0 0 0 0 — 


(b) One Wall Straight 


The equations appropriate to this case are (29), (23) and (30), the first of which 
is satisfied identically in virtue of the symmetry about ¢=0. Equation (23) is used, 
as in (a), to define B,. The assumption HU=1-68396a (selected in an effort to 
obtain R/H nearly equal to 1/2), and the division of the range into sixteen equal 
parts, leads eventually to the solution set out in Table V. For this case it is found 
that R/H=0-4749. The last entry in column 5 is not reliable as it changes rapidly 
with change in 2; in fact increasing z from 19-81° to 20-31° reduces ¢ from 0-227 to 
0-198. Some small error in « at ¢/(4a)= 15/32 is an inevitable consequence of using 
the approximation (27) to calculate sU /(4a) at this point. The same remarks apply 
to the last entry in column 5 of Table IV. 


The standard result for the first order term in solid blockage is 


~ 


which may be established from equation (3) by an argument parallel to that already 
given for constant pressure walls. For the circular cylinder 


—0-0927, 


which is only about 37 per cent. of the average result obtained from column 5—in 
marked contrast to the result obtained for constant pressure walls. Fig. 5 shows 
q/U and « as functions of <. 


The results obtained for the blockage factor lead to the conclusion that, as 
far as drag corrections are concerned, relatively larger models can be safely used 
in open-jet tunnels than in closed tunnels. 
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Fig. 5. 
Solid blockage and velocity distributions for a circular cylinder. 


10. The Setting of “Streamline” Walls 


Equation (8) has an important application to the theory of flexible-wall wind 
tunnels. In such tunnels the wall shape can be changed as required, and the 
principle adopted is to set them, as nearly as possible, to lie along streamlines of 
the infinitely deep stream. This eliminates any need to apply wind tunnel 
corrections, and, more important, avoids that defect of straight-walled tunnels at 
high subsonic Mach numbers known as “ choking.” 


The method of setting these walls is shown in Fig. 6. In this figure ABC is the 
constant pressure wall in the absence of the aerofoil, A,B,C, is the constant 
pressure wall in the presence of the aerofoil (assumed to be at zero incidence on 
the tunnel centre line), while A,B,C, is the corresponding free streamline about the 
aerofoil. 
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The setting of streamline walls. 


These curves have a common point at an infinite distance upstream, but 
downstream the existence of the aerofoil’s wake causes a displacement of the 
constant pressure and free streamline walls from the wall position of the empty 
tunnel. The constant pressure wall is obtained experimentally, so the displace- 
ment BB, at any point will be known. The required streamline wall A,B,C; can 
now be determined from the equation BB,=«BB,, provided that the value of « is 
known. The present purpose is to determine the theoretical value of this ratio «. 


Present practice, based on calculations of Lock and Beavan", is to take « as 
constant and equal to 0-6. The calculations of Lock and Beavan (for doublets, 
vortices and sources instead of aerofoils) are not consistent with the method 
actually used in practice for setting streamline walls. In their calculations they 
made the constant pressure and streamline walls pass through the same point 
opposite the doublet (source or vortex), whereas in experimental practice the walls 
pass through a common point effectively an infinite distance upstream (x= — 00). 
Thus for the doublet their method results in values of « of 4x and 16/=* (+0516), at 
x= -—00 and x=0 respectively, whereas if the walls have a common point at an 
infinite distance upstream it is shown later that the values of « are 00 and 2/= 
(=0-637) respectively. Furthermore it is shown later that in the presence of a wake 
*=1 at an infinite distance downstream, whereas Lock and Beavan recommended a 
factor of 0:6 throughout the range of x. The effect of this is not serious when 
applied upstream, but downstream, where the wake produces a permanent displace- 
ment of the tunnel wall from the empty tunnel position, to return the wall 0-4 of 
the displacement will plainly result in the reappearance of some of the wake 
blockage. This may not, of course, materially affect the experimental results. 


It is suggested in this paper that for symmetrical aerofoils at zero incidence « 
should be taken as 2/= along the wall opposite the aerofoil, and then parabolically 
increased such that it becomes unity at 38. tunnel widths up and down stream from 
the centre of area of the aerofoil. Outside this range the constant pressure walls 
could probably be left undisturbed (« =1) with little resulting error. The analysis 
leading to these conclusions is as follows. 


By considering, as in Section 9, only that half of the flow for which Y >0, the 
flow can be regarded as being one in a “channel,” the aerofoil surface being y=0, 
and the constant pressure or streamline wall being Y=h/m,. The equations 
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appropriate to this case are (5), in which H is now half the tunnel height, and (8). 
In the limit, H —> 00, these equations yield 


which applies to a symmetrical aerofoil in an infinite stream. On the streamline 
y=h/m,, the imaginary part of this equation gives 


00 


which enables the shape of the streamline wall to be calculated directly. The 
approximation of thin aerofoil theory is now made by replacing 9,, 9’ and ¢ in (34) 
by dy,/dx, Ux’ and Ux. Then from (5) and (34) 


2’ = 


where c is the chord length, and », is the displacement of the wall Y=h/m,, 
from its original straight position in the empty tunnel. The origin of x is taken at 
the mid-chord point, and the range of integration is reduced, as ,=0 outside 
Now y,;=0 at x= so integrating (35) gives 


—c/2 


Similarly it is found that on Y=h/m,, equations (5) and (8) yield the displacement 


a’ = —c/2 


for the constant pressure wall (cf. the third equation after (27)). 


Now y,=0 at x=-—4c (the leading edge), but at the trailing edge (x=4c) 
Yo=}cCp, where Cp is the drag coefficient. Here it has been assumed that a 
“sting” of thickness 4cCp is attached to the trailing edge to represent the wake 
displacement thickness (see the theory in Ref. 3). 
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For large z, 


gd (2)= z— sech*z..., and tan- 


and it therefore follows from (36) and (37) that well downstream (x=00) both y, 
and y, tend to }cCp, that is «(=y;/y,) tends to unity, whereas well upstream o 
tends to infinity like (cosh z)/z. 


The value of ~ is now calculated in the neighbourhood of the aerofoil. From 
(36) and (37), 


2’ = —c/2 


2 
For small z, 


gd (2)= 52 (1- tan-? (2)=2(1- 


so that for small values of x/H we can write 


2’ = —c/2 


Integrating by parts and ignoring terms of order (c/H)’, it is found that 


-1) 


where J (x) is the second moment of the profile section about x. It is clear from 
(38) that.o will always exceed 2/x, but not by very much in the neighbourhood of 
the aerofoil. Whitehead‘ has also obtained the first order result c=2/ by apply- 
ing conformal transformations to special aerofoil profiles. 


It is of some interest to determine the value of x at which « becomes unity. 
From (36) and (37) this occurs where 


e/2 


[r(2 


2’=—-—c/2 


dy, (x)=0, . 689) 


where F (y)=gd(z/2)y—tan-'y. This equation can be solved numerically for a 
given aerofoil if Cp is known. A good estimate of the value of x can be found 
as follows. . 
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For the rectangle y,=43t, —4c<x<4c, relation (39) gives 


F (244) - (1-262) ‘(4% (40) 


which is easily solved graphically. When cC,/(4t) can be ignored, corresponding 
values of-c/(H8,,) and x/(H8,,) satisfying (40) are 


c/(HB,.) 0 0:2 0-4 0-6 0:8 1:0 1-2 
x/(HB,.) +143 4+1:44 41:45 +146 +1:48 41:51 41°54 


This shows that x/(H8,,) is practically independent of c/(H§,,) and, since an aero- 
foil can be constructed from a number of rectangular strips of varying chord, it is 
reasonable to apply this table to any aerofoil shape. In fact even for the extreme 
case of the cylinder of Section 9.2, where c/(H8,,)=2 x 0:5278, the table yields 
x/(HB,,)= 1-52, a value less than one per cent. in error from the true value of 1-532 
(see Fig. 4). If Cp is not made zero the values of x/(H8,.) are reduced downstream 
and increased upstream, but the effect appears to be small for typical values of the 
ratios c/(H8,.) and cC)/(4t). Thus it would probably be a satisfactory rule in 
experimental work to take o as unity outside the range —}8,,<x/(2H)<?f,. 
say, as the adjusting of the constant pressure wall to the “streamline” position has 
to stop somewhere upstream and downstream, and the points where »=1 are 
convenient for this purpose. 
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Effect of Yielding on the Fatigue Properties 


i | of Test Pieces Containing Stress 

Concentrations 
ling f 

K. GUNN, B.Sc., A.R.S.M. 

(Aluminium Laboratories, Banbury) 
—— i Summary: A method of estimating the fatigue properties of notched test 

pieces when cyclic loads are applied in the presence of mean stresses is 

-rO- | described. The method takes into account the local yielding which occurs 
itis | near the base of the notch when the nominal stresses are sufficiently high. 
sme ) The general form of the predicted R/M diagrams closely resembles those 
F id obtained experimentally in previously published work. 
elds | 
532 
2am 1. Introduction 
the | A large number of investigations have been made in the past on the effect of 
> in 


mean stresses on the fatigue properties of aluminium alloys. The published reports 
Box. ) of these investigations have not formulated general relationships, since for the most 

| part they have been merely a description of test results. An attempt is made here 
to rationalise the results of tests on notched test pieces which are subjected to mean 
stress by postulating simple relationships between the nominal applied stress and the 
actual stress present in the test pieces. 


; A description of the derivation of a general relationship applied to idealised 
sion | materials is given and the predicted values of fatigue strength are compared with 


ibed | actual results obtained by experiment in work previously published by other 
onal investigators. 
the 
NOTATION 
R__ semi-range of stress under mean stress M 
and : R, — semi-range of stress under zero mean stress 
M mean stress 
eigh- ff U___ ultimate tensile stress 
K, elastic stress concentration factor 
The K, “plastic” stress concentration factor 
E  Young’s modulus 
om E, secant modulus at base of notch 
0. Received June 1955. 
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E,, secant modulus at a point away from notch but in weakest section of 
test piece 


j fatigue strength of plain test pieces ) 
h reduction factor\ =~ 
fatigue strength of notched test pieces 
_ K,-1 

Y _ yield stress 

f, | nominal stress in weakest section of test piece 

f. local stress at base of notch 


€,, €, total strains corresponding to f, and f, respectively 


2. Initial Assumptions Based on Modified Goodman Relationship 


For simplicity it may be assumed that when plain test pieces are subjected to 
fatigue loads in the presence of mean stresses, the relationship between Range of 
Stress and Mean Stress for a given endurance is 


R=R(1- ‘ ‘ (1) 


This relationship is known as the modified Goodman relationship. (It is not 
necessary to assume that this relationship is correct; it does, however, make the 
development of the theory more easily understood). 


In a notched test piece under static stresses which are entirely elastic, the 
maximum stress at the base of the notch (termed “local” stress in this paper) is 
given by 


fe=f, x Ky. 


It is assumed for the present that this stress concentration factor K, can also 
be applied to fatigue test results. The nominal semi-range of stress under zero 
mean stress which can be applied to the test piece for a given endurance will be 
reduced from R, for plain test pieces to R,/K;, for notched test pieces. The mean 
stress at the base of the notch will be increased over the mean stress M in the body 
of the test piece by the factor K, and will equal K.M. These values, substituted in 
equation (1), give for notched test pieces the equivalent relationship to the 
Goodman rule :— 


The equations (1) and (2) are illustrated diagrammatically in Fig. 1. This diagram 
is here called an R/M diagram. 
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Ry 
« 
PLAIN TEST PIECES 
es} 
| R NOTCHED TEST PIECES 
) 7% ne Be 
z 
° 
i z 
u 
hip MEAN STRESS = M 
d to Fig. 1. 
eof | Modified Goodman relationship for plain and notched 
test pieces in the absence of local yielding. 
(1) | Equation (2) is only applicable when a test piece is stressed purely elastically, 
) for once yielding occurs the relationship f,=f, x K, is no longer valid. The effect 
on equation (2) of local yielding in notched test pieces is discussed in the next 
h section, where it will be seen that higher values for the nominal range of stress than 
and those predicted by equation (2) can be tolerated, once yielding occurs. 
the 
vis | 3. Effect of Yielding on Mean Stress 
! In a notched test piece under cyclically varying loads, the material which yields 
. at the base of the notch will be hardened sufficiently to prevent further plastic 
deformation* (see Appendix I). (The hardening may be completed in the first cycle 
or in several cycles; the actual number of cycles will depend on the properties of the 
also test piece and the applied stresses). When the test piece has “ settled down ” and is 
zero stressed purely elastically the stress amplitude at the base of the notch will have 
| be | the same value as it would have had if the body had been purely elastic from the 
lean beginning (see Appendix I). Although the actual range of stress at the base of the 
ody notch will be equal to K, times the nominal range of stress, the local maximum 
d in stress will be less than the stress derived from purely elastic considerations because 
the the local maximum stress will be limited by yielding, its actual value depending on 
} the shape of the stress-strain curve. As a result of the diminution of the local 
; maximum stress, but the maintenance in full of the local range of stress, the 
a | local mean stress in the cycle will be reduced below the value 
K, x nominal mean stress. 
*The material is assumed to be a homogeneous isotropic solid and plastic deformation under 
apparently purely elastic stresses is not considered. 
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With this relief of local mean stress the range of stress which can be withstood under 
a given nominal mean stress will be greater than that predicted if local yielding had 
not been taken into account. 


4. The R/M Diagram of Notched Test Pieces 


Using the concept that due to local yielding the local mean stress diminishes, it 
is possible to derive R/M diagrams for notched test pieces. The method based on 
the theory outlined in the previous section is presented in this section. 


Under purely elastic stresses the range of stress, related to the mean stress by 
the equation R=(R,/K.)[1—K.(M/U)], decreases with increasing mean stress 
along the line AB in Fig. 2. After point B the maximum stress in the local 
stress cycle at the base of the notch exceeds the yield point of the material and so 
the local stress is no longer directly proportional to the nominal stress in the body 
of the test piece. 


The stress concentration factor after yielding has occurred can be calculated 
(see Appendix II) from the following formula 


Using values of K, derived by means of this expression, the local maximum stress 
can be calculated from each chosen nominal maximum stress. From the now known 
local maximum stress the corresponding local range of stress can be obtained from 
Fig. 3, which relates the semi-range of stress and the maximum stress. This diagram 


PLAIN TEST PIECES 


NOTCHED TEST PIECES 


eNOMINAL SEMI-RANGE OF STRESS~-R tons/sq. in. 


NOMINAL MEAN STRESS - M_ tons/sq- in 
Fig. 2. 


The influence of local yielding in notched 
test pieces on the R/M diagram. 
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is easily derived from the modified Goodman relationship R=R,(1—M/U) 
illustrated in Fig. 1. The nominal semi-range of stress and the nominal mean stress 
in the body of the test piece are obtained directly by 
Local semi-range of stress 

K, 


Nominal semi-range of stress= 


and 


Local semi-range of stress ) 
K, 


Instead of using the “semi-range/maximum stress” diagram. of Fig. 3, the 
local range of stress can be obtained from the diagram in Fig. 4, where the maximum 
and minimum stresses are plotted directly against the mean stress. The nominal 
range of stress and the nominal mean stress can be derived as before. 


Nominal mean stress=Nominal maximum stress -( 


The post-yielding relationship between nominal mean stress and nominal 
semi-range of stress (calculated in this manner), is illustrated in Fig. 2 by the curve 
BCD and in Fig. 4 by BCD and B’C’D. The curve BC in Fig. 2 (or the curves 
BC and BC’ in Fig. 4) illustrate the relationship when the body of the material is 
elastic and only the material at the base of the notch undergoes plastic deformation. 
At C in Fig. 2 (or C and C’ in Fig. 4) the body of the test piece also begins to yield 
and this causes the corresponding range of stress to diminish rapidly as the ultimate 
stress is approached. 


5. Secondary Modifications to the Theory 


It is obvious that the idealised curves derived from the simplest assumptions 
can be modified to cover a number of additional factors. Some of these factors are 
now discussed. 
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Fig. 4. 
The influence of local yielding on the fatigue properties of notched test pieces. 


5.1. MODIFIED GOODMAN RELATIONSHIP 


The initial assumption which was mentioned earlier was only included to 
demonstrate the theory more easily. Any relationship between R and M for plain 
test pieces can be used, since the shape of the R/M diagram for notched test pieces 
is determined by graphical methods. 


5.2. STRESS CONCENTRATION FACTOR 


The expression in equation (3) for the “ plastic” stress concentration factor K, 
cannot be solved directly for any chosen nominal stress. Either a graphical method 
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Stress-strain curve used in deriving Fig. 2. 


or a process of successive approximations must be used once yielding has occurred. 
A simpler, but less accurate, conception is to assume that the strain concentration 
factor is equal to the elastic stress concentration factor after local yielding, as well 
as before local yielding occurs. A value for the local maximum stress can be 
obtained directly from the stress-strain curve as shown in the example in Fig. 5. 
In this example a nominal stress of 17 tons/in.* was chosen—its associated elastic 
strain being 0-36 per cent. With an elastic stress concentration factor of K,=3-0, 
the local strain is taken as 3 x 0-36 per cent.=1-08 per cent. The stress correspond- 
ing to this strain is 21:2 tons/in..*. The semi-range of stress, 4:35 tons/in.’, 
associated with this maximum stress can be read from Fig. 3. The nominal 
semi-range of stress is 1-45 tons/in.* and the nominal mean stress is therefore 
15-55 tons/in.?. These values can be compared with the curve in Fig. 2, where it 
can be seen that the error due to the simpler method is small. 


A comparison of the stress concentration factors calculated by both methods, 
using the stress-strain curve in Fig. 5 and an elastic stress concentration factor 
K,=3-0, is shown in Fig. 6. Only at high stresses is there appreciable divergence 
in the two curves and so the simpler method has been adopted in the examples 
illustrating the theory. 


5.3. STRESS AT WHICH GENERAL YIELDING BEGINS 


The shape of the curve CD in Fig. 2 is not accurately defined. The first 
uncertain part is the point at which general yielding occurs. Often with notches 
there are triaxial stresses induced in the material away from the notch, as for 
example at the centre of a cylindrical test piece containing a circumferential notch, 
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Variation of stress concentration factor when local 
stresses are greater than the yield point. 


Two expressions for K, are illustrated. 


Under triaxial stresses, assuming all are in tension, the nominal yield stress will rise 
and therefore there will be a prolongation of the flattened portion BC of the 
R/M diagram. This may be partly an explanation of the large range of stress 
which is sometimes obtained experimentally under extremely high mean stress as in 
Fig. 9(b) (this figure is discussed later). 


5.4. THE R/M DIAGRAM UNDER HIGH MEAN STRESSES 


The high value for the points under high mean stress discussed in Section 5.3 
is but one example of the uncertainty of this part of the diagram. The precise 
shape of the R/M diagram for plain test pieces under very high mean stresses is 
virtually unknown. Indeed it is likely to depend on the testing technique employed 
in the fatigue tests. The application of static test results or dynamic test results 
from plain test pieces under mean stresses approaching the ultimate stress, to a 
prediction of the behaviour of notched test pieces in fatigue is therefore unreliable. 


It is possible that true stress-strain curves derived from tests on notched test 
pi¢ces, together with the true stress-strain curves obtained on plain test pieces, could 
be used in determining the maximum local stress in a loading cycle, but no attempt 
has been made to do so in this paper. 
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Fig. 7. 


(d) K,=3°4. 10’ cycles. 


Comparison of theoretical R/M diagrams for notched test pieces and experimental results. 
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Comparison of theoretical R/M diagrams for notched test 
pieces and experimental results. 24S-T4®), 


6. Experimental Verification of the Theory 


Test results taken from Lazan and Blatherwick™, and Wallgren’, are illustrated 
in Figs. 7 and 8 in the form of R/M diagrams. Fig. 7 shows the results of fatigue 
tests on both plain and notched test pieces in 14S-T6 (an alloy similar to those con-  _[ 
forming to B.S. L65) subjected to various mean stresses. Two designs of notched | 
test pieces were tested and the fatigue strengths for two endurances, 10* and 
10’ cycles, are given for each notch. Figs. 7(a) and 7(c) illustrate the results for a 
notched test piece with an elastic stress concentration factor of K,=2-4 and 
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(b) 
Fig. 9. 
Comparison of experimental test results and theoretical curves assuming 
(a) K,=3-4 and (b) K,=2.6. 
Alcoa 24S-T4. 10? cycles endurance”). 


Figs. 7(b) and 7(d) illustrate the results for a notched test piece with an elastic 
stress concentration factor of K.=3-4. Similarly in Figs. 8(a) and 8(d) the results of 
plain and notched (K,=2-05) test pieces in .24S-T4 are illustrated for 10* and 
10’ cycles endurance. Superimposed on the diagrams are calculated curves derived 
from the results from plain test pieces and the elastic stress concentration factors. 


A striking feature of all the experimental results is the presence of a more or 
less well defined knee similar to that shown by the calculated curves. In most 
cases however the calculated curves fall below the experimental results, because in 
deriving the curves it was assumed that the elastic stress concentration factors would 
exert their effects to the full. It is a well known experimental observation that this 
is generally not so and that the strength reduction factor K, is often less than the 
elastic stress concentration factor K;. A relation between these two factors has 
been expressed“ as the “sensitivity index” q defined by 


K,-1° 
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The sensitivity index is often called the notch sensitivity factor. The numerical 
value of q depends, for any alloy, on the condition of heat treatment, on the test 
piece size and shape, on the magnitude of the mean stress and on the number 
of cycles endurance. By ignoring g in the development of the theory it has been 
tacitly assumed that its numerical value is always unity. This can cause the 
effective stress concentration factor to be overestimated (by unknown amounts) with 
consequent discrepancy between the experimental and calculated curves. 


If instead of K, the value K, calculated from the tests results under zero mean 
stress were used, it would for example modify the equation (2) so that it became 


It would also modify the R/M diagrams for notched test pieces obtained graphically; 
for example Fig. 9(a) shows experimental results’) (plotted in the form of R/M 
diagrams) for plain and notched test pieces in 24S-T4, together with the calculated 
R/M diagram obtained by using K;, to its full extent, i.e. K,=3-4. In Fig. 9(b) how- 
ever, the calculated R/M diagram obtained by using the value K,=2-6 instead of 
K, is compared with the same experimental results. In this figure the experimental 
results are seen to lie close to the theoretical curve, especially at low and inter- 
mediate mean stresses. 


7. Simplified Procedure 


The derivation of the complete R/M diagram for notched test pieces entirely by 
graphical methods is both complicated and tedious. A simplified procedure is 
outlined here which is rapid when graphical methods are used and is amenable 
to algebraic treatment. 


First, the R/M diagram for plain test pieces is drawn. The first part (i.e. before 
local yielding) of the R/M diagram for notched test pieces (AB in Figs. 2 and 10) is 
drawn from the plain test piece curve by reducing both the semi-range of stress and 
the mean stress by the elastic stress concentration factor K;. The point of local 
yielding is obtained by the intersection of this R/M curve and the straight line 
joining the points “stress=(yield)/K,” plotted on both axes (i.e. the straight line 
that is the locus of maximum stress=(yield)/K;,). “ Yield,” itself, is taken as the 
0-1 per cent. proof stress in aluminium alloys. 


It is now assumed, to simplify the diagram, that the range of stress does not 
change with increase in mean stress after local yielding has occurred (line BC in 
Fig. 10). There is little error in this assumption until general yielding is approached, 
but as such high stresses are outside the scope of normal design it is convenient to 
ignore this region of the diagram. The simplified diagram is shown in Fig. 10 and 
under a different method of plotting in Fig. 11. The R/M diagram obtained by 
considering work-hardening (i.e. the curve BC in Fig. 2) is superimposed on these 
figures to illustrate the small degree of error which the simplification entails. 
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SIMPLIFIED R/M DIAGRAM 


-—<—-— CALCULATED R/M DIAGRAM 


YPLAIN TEST PIECES 


NOMINAL SEMI-RANGE OF STRESS — tons/sq in 


20 30 
NOMINAL MEAN STRESS — tons/sq.in. 


Fig. 10. 
A simplified R/M diagram. 


R 4/3 
AB: (x, 3) 
K, U 


R E 


BC: R=— 
K, 

The algebraic treatment is straightforward, since the R/M diagrams for plain 

test pieces can be expressed by a number of equations, one group being of the form 


R=R,[1-(#)']. 


Equation (1) is of this form, n being taken as 1. Similarly equation (2) can be 
obtained by putting n=1 in equation (6), which is the equation for the first part of 
the R/M diagram for notched test pieces : — 


The intersection of this curve and the straight line representing 


maximum stress=(yield)/K, 


Y 
ie. M+R ‘ ‘ (7) 


gives the start of local yielding. This point (and hence the horizontal line BC in 
Fig. 10) is obtained by substituting equation (7) in equation (6). 
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CALCULATED R/M DIAGRAM 
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Fig. 11. 
Simplified maximum and minimum stress diagram. 


R 4/3 
AB and A’B’: R= | 
K, U 


BC and B’C’: R= 
K, 


Y-K,R 


290 


30 


The Aeronautical Quarterly 


| 
© 

= / 

7 

Je 
ite) 
— | 
® 
® | 
! 
10 
= 


irterly 


FATIGUE OF NOTCHED TEST PIECES 


The curves for these equations, with the constant n taken as 4/3, are illustrated in 
Figs. 10 and 11. 


In the previous section an example was given of the use of the strength 
reduction factor K, in the place of the elastic stress concentration factor K, A 
word of warning is necessary against the unreserved use of this procedure in practice. 
There is evidence to suggest that the notch sensitivity factor g approaches unity 
(i.e. the strength reduction factor K, approaches the elastic stress concentration 
factor K,) as the volume of highly stressed material increases. If the factor K, were 
obtained from test pieces of the same order of size as a service member, then it is 
possible that it could be used with safety. When the service member is large com- 
pared with the test pieces, then it is advisable to use the elastic stress concentration 
factor to its full extent. A further warning is required against an optimistic inter- 
pretation of plain test piece results, since the location of the point of local yielding, 
occurring as it does at the intersection of two steeply sloping lines, may be critically 
dependent on minor changes in shape of the initial R/M diagram. 


This application of K, or K, disregards other factors which may be important 
in promoting the onset of fatigue failure. For example, if corrosion or fretting 
corrosion cannot be ignored then they must be taken into account by multiplying 
the stress concentration factor K, by an empirical factor. Apart from these adverse 
service conditions it is expected that it may be difficult to determine the value for 
the stress concentration factor itself; however, it is necessary to have specific details 
of designs before this can be discussed. 


8. Discussion 


In the development of the theory it has been assumed that the static stress- 
strain curve adequately describes the behaviour of the material under dynamic 
loads. Insufficient information is available at present to enable a more accurate 
estimation of the stress-strain relationship under fatigue conditions, but it must be 
borne in mind that a dynamic stress-strain relationship may modify the final R/M 
diagrams appreciably. 


However the qualitative agreement between the theory and actual test results 
justifies further speculation on the behaviour of notched test pieces. It is evident 
that the theory lends some support to the belief that for good fatigue properties a 
large difference is needed between the proof stress and the ultimate tensile stress (as 
in material where the ratio (proof stress)/(ultimate tensile stress) is low). The 
redistribution of stress near the base of a notch will be more easily accomplished 
in such a material than in a material with a high ratio (proof stress) /(ultimate tensile 
stress), and so under certain mean stresses it is possible that the former will have 
relatively higher fatigue properties than the latter, purely because of the difference in 
yield strength. It is necessary, however, that this suggestion be supported by 
experimental evidence before it can be accepted. 
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Appendix I 


A test piece containing a notch, when subjected to tensile stresses, will eventually 
undergo plastic deformation round the notch, as in the shaded area in Fig. 12. The 
local stress at the base of the notch will be given by the point A on the stress-strain curve 
of Fig. 13, when the nominal stress is increased from zero (point O) to point P. On 
releasing the load, if the volume of material plastically deformed is very small compared 
with that still purely elastic, the local stress will be B and the stress in the body of the 
test piece O. On similar subsequent applications of tensile stress the local stress at the 
notch will cycle between B and A, i.e. purely elastically but with a lower mean stress 
than that in the body of the test pieces multiplied by the stress concentration factor. 


It is possible that yielding in compression could occur round the notch on releasing 
the load. In this case, several cycles must elapse before the local stresses are purely 
elastic, but the mean stress is still less than that given by K, x the nominal mean stress. 

With zero nominal mean stress the local mean stress will also be zero if the 
properties in compression are the same as in tension. 


O 


(b 
(Q) FIRST HALF CYCLE (b) SECOND HALF CYCLE} 
Notched test pieces under tensile stress. Stress in body of test piece: P and O. 
(a) Purely elastic. (b) Partially plastic Local stress round notch: A and B. 


(K, x nominal stress is greater than the yield 


point). 
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Graphical solution of the equation f,/f,=1—(K,—1) e,/e, when K,=3-0. 


Appendix II 


The equation extracted from Ref. 1 relating purely elastic stress concentration 
factors with the partially plastic stress concentration factor is 


E, 
K,=1+(K,- 1) (3) 
oo 
E, 
This can be re-written — =1+(K,- 1)}—. (9) 
E,. 
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This expression can be solved in three stages : — 


(i) When the material at the base of the notch is still elastic, the elastic stress 
concentration factor K, operates to its full extent and so curve AB in Fig. 2 is valid. 


(ii) When the material at the base of the notch has yielded locally but the body 
of the material is elastic, equation (9) becomes 


f E, 
7, =1+(K,- 


By taking values for f, and reading E, from the stress-strain curve as in Fig. 5, f, } 
can be determined in terms of K, and f,/f, can be calculated. 


(iii) When the whole material has yielded, a graphical solution is necessary since 
the relationship between stress and strain is not expressed algebraically. The equation } 
(3) between K, and K, can be re-written 


f, I, ¢ 
7=14+(K,-1)- +. 
‘ 
This reduces to —=1-(K,-1)—. 


By taking values for f, and the corresponding values for e, from the stress-strain curve 
of Fig. 3, both sides of the expression can be plotted against values of e, as in Fig. 14. 
From these curves we can obtain values of f,/f, (i.e. K,) which satisfy the equation. 


The values so obtained for K, are shown in Fig. 6. 
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: Stabilisation of the Aeroplane in Symmetric 
dy Flight at Zero or Very Small Air Speeds 


W. S. COLEMAN, B.Sc., Ph.D. 
(Blackburn & General Aircraft) 


we 


Summary: The possible use of jet thrust as a source of lift, either by 
deflection of the force from the normal propulsion axis, or by some other 


wd means, raises a new stability problem under conditions of hovering 


on flight and flight at speeds for which the aerodynamic damping becomes 


insignificant. At the same time, despite the small dynamic pressure, the 
direct air loads cannot be neglected, for, due to the very large positive or 
negative incidences peculiar to the motion, the associated coefficients are of 
a compensatingly high order. Consequently, the problem retains an aero- 
dynamic element, but is principally one of automatic stabilisation, amenable 
to treatment by ordinary mechanical or electrical circuit theory. The investi- 
gation has been restricted to the longitudinal-symmetric motion with virtually 
the only workable form of automatic control, namely, one responsive to the 
angular displacement and velocity in pitch, but these conditions undoubtedly 
embrace all the essentials of any practical application, e.g. in take-off and 
landing by the new technique. The solution to the stability in a gust, or ina 
succession of disturbances, may then be easily obtained in closed analytical 
form, but the more general case of unsteady motion, when the displacements 
i are no longer purely oscillatory, necessitates, in part, the use of step-by-step 
integration. This arises from the added complexity of the aerodynamic terms, 
due to the appearance of significant changes of incidence, and to the very 
large and variable values of the resistance derivatives at the angles involved, 
so that, in the present state of knowledge, the forces determining the motion 
cannot be expressed functionally. The labour which such calculations entail 
may be considerably shortened, however, by the application of a procedure 
based on mean deviations, but it is also found that the conclusions to be 
drawn from the hovering case remain valid in a number of important respects 
for the general problem. This feature is demonstrated in the discussion at the 
end of the paper on a number of examples designed to illustrate the charac- 
teristics of the automatic control and the conditions necessary for stability 
in a variety of circumstances, both when the aircraft is hovering and when 
its displacement increases continuously with time. 


rve 
14. 
on. 


1. Introduction 


The sustentation of an aeroplane in the normal flying attitude by means of the 

thrust from a turbo-jet engine, or a number of such engines, is a development which 

*» will probably materialise in a comparatively short time. To obtain the necessary 

|} vertical reaction to enable the aeroplane to hover, and to rise and fall, either the 

engines intended primarily for normal propulsion may be used, in which case the 

|| jets must be temporarily deflected, or a battery of small power units, installed 
exclusively for the production of the necessary lift, might be adopted. 
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As a result of the possibilities which this development offers, much thought is 
being devoted to the problems involved, both in relation to the engine, and to the 
control and performance of the aeroplane, under the new conditions which are 
introduced. The present paper discusses the stability during the stages when the 
mean air speed is so low that the aerodynamic damping of a disturbance is negligible. 
Accordingly, stabilisation of a mechanical character is essential, and it is assumed 
that when the aeroplane is airborne in a condition for which the aerodynamic control 
and stability are ineffective, there is no intervention on the part of the pilot, the 
motion being determined purely automatically in terms of the displacements and 
velocities. Again, for the purpose of the investigation, it is assumed that the system 
concerned consists of the necessary servo-mechanisms coupled to subsidiary jets 
disposed at suitable points on the aeroplane. In considering the stability under these 
conditions, it is also clearly sufficient to stipulate the form of the response without 
further examination of the detailed functioning of the auto-controls, and in this 
connection it is assumed that there is a linear correspondence between the input 
and output signals. At the same time, lag in the system as a whole is taken into 
account. Lastly, it is manifestly out of the question, within the space of a short 
paper, to consider more than one or two of the very wide range of ways in which 
such controls may be designed to function through various measures of the 
disturbances to which the aircraft is subject. Accordingly, attention is confined here 
to what would seem to be the two most practical, and commonly used, components, 
namely, the angular displacement and the angular velocity which, in the present 
instance, refer essentially to the pitching motion. 


Although, as discussed in Appendix I, all the relevant aerodynamic damping 
terms may be ignored in the immediate problem, this is not the case with respect 
to the corresponding static forces and moments for, despite the very small air speeds, 
the incidence is clearly large, tending to the limit +90°, and wind tunnel measure- 
ments of the lift, drag and pitching moment on a model of a typical high performance 
aeroplane under such conditions indicate that the coefficients then become so large 
as to make the actual forces and moments significant. Indeed, if the degree of 
stability is inadequate, the aircraft rapidly executes a loop into the ground, due to 
the combination of the thrust and the aerodynamic terms (see Section 4.4). 


It will also be assumed that the motion is confined to the longitudinal-symmetric 
state, since this clearly accords in all essentials with the relevant and normal flight 
conditions. It implies, however, that any gyroscopic yawing couple generated by 
the rotating parts of the power units and other like sources, when the aircraft pitches, 
is either negligible or else so damped as to produce no significant departure from 
two-dimensional motion. This point is considered in detail during the early stages 
of a disturbance in Appendix II, from which it appears that, for any reasonable 
magnitude of the couple, the deviations from the truly longitudinal-symmetric state, 
provided there is adequate stabilisation in both pitch and yaw, will be sufficiently 
small to satisfy the foregoing assumption. 


Because of the small linear velocities which are dealt with, both as regards the 
horizontal and vertical components, the problem of the stabilisation associated with 
vertical flight may be largely studied in terms of the hovering aeroplane subject to 
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a gust of known characteristics. It will be seen that all the essentials of the motion 
during a typical take-off on jet lift, for example, are exhibited under these simple 
conditions which yield a solution of the equations of motion in closed analytical 
form. The results also isolate the effects of the relevant geometrical and inertial 
characteristics of an aircraft in a very simple manner. On the other hand, if the 
aircraft is not hovering, and the resultant mean air speed is, in general, an unknown 
function of the time, the motion can only be determined from a step-by-step solution 
for each individual case. The labour inherent in such calculations has been much 
reduced, however, by the mean deviation method fully described in Section 3.2. 


The first part of the paper, therefore, is devoted to the hovering case with a 
variety of disturbances, and these results are amplified later by a selection of the 
more interesting conclusions from a number of step-by-step analyses of the general 
case when the displacements of the aircraft are not purely oscillatory. 


There remains one other point to consider at this stage. Because the flight 
path of the aircraft is a primary consideration in the present investigation, fixed 
axes of reference are more convenient than the customary moving axes of stability 
theory. At the same time, owing to the symmetric character of the motion, freedom 
in rotation is confined to pitch, and hence the corresponding inertial term requires 
no resolution, since the transverse body axis and its counterpart in the fixed 
reference system remain normal to the plane of symmetry, and therefore parallel to 
one another. Accordingly, the stationary reference frame has been adopted and the 
appropriate forms of the equations of motion are considered with respect to it. 


NOTATION 


x, z fixed rectangular co-ordinates*, with the x-axis horizontal, the 
z-axis vertical and the origin at sea level 


t time 
7 dummy time variable 
T _ period of oscillation of the aircraft 
W,m_ weight and mass of the aircraft respectively 


B- moment of inertia of the aircraft with respect to the transverse 
axis (normal to the xz-plane) through the centre of gravity 


P total thrust developed by the aircraft 


M total applied pitching moment with respect to the centre of 
gravity 


L, D, M aerodynamic lift, drag and pitching moment about the centre of 
gravity respectively 


Ci, Cp lift and drag coefficients 
C,, moment coefficient 
o relative air density 


*It is convenient to regard the reference axes as fixed relative to the earth, and the motion, 
apart from the gust conditions investigated, as occurring in still air. 
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B 
Pp 
lw z 
Fig. 1. 
Co-ordinate system. 
fp, Standard sea level air density 
V true air speed 
S gross wing area 
C geometrical mean wing chord 
a angle of incidence with respect to either wing chord or longi- 
tudinal body axis, the two being assumed parallel in the present 
investigation 
B _ slope of the flight path at any point with respect to the horizon- 
tal, taken positively in a climb 
y inclination of the jet pipe axis to the longitudinal body axis, 
assumed positive when measured downwards from the zero 
position for normal propulsion 
6 angular displacement in pitch relative to the horizontal, 
i.e. 
w vertical gust velocity 
Ww, maximum value of w 
T, duration of a gust 
o = t/T, 
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stiffness coefficient of an automatic control which responds 
linearly with 6, i.e. the moment generated per unit of angular 
velocity in pitch 

stiffness coefficient of an automatic control which responds 


linearly with 9, i.e. the moment generated per unit of angular 
displacement in pitch 


= —K,/(2B) 


/(—K,/B) 

(s?—a’) 

(2B) 

«,/(1—2ae,) 

time lag of a 6-control system 
time lag of a 6-control system 
(2a —s*e,)/ {2 (1 —2ae,)} 
s//(1—2ae,) 


= (s,?—a,”) 


.. 
= 


C, o C, 
C, 


2, 
Fit) 
hy, hy, hy, hs 


Bo» 81> 82> Bs 
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Laplace transform { be-*dt 


unit step function 

AM/B 

amplitude of a square wave pitching-moment disturbance 
half excitation period of such a disturbance 

a? — b? + 4? 

a’ + b? —4w? 

a,? —b,?+ 40? 

a,2+b,?—40? 

constants of integration defined by equation (17) 
constants of integration defined by equation (24) 


constants of integration defined by equation (17) when the fol- 
lowing substitutions are made, namely, C,’ for C,, C,’ for C,, 
etc., and a, for a, b, for b, x, for x,, A, for A, v, for v 


constants of integration defined by equation (36) 
the variable M/B 
coefficients in the step-by-step expansion of x 


coefficients in the step-by-step expansion of z 
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2. The Equations of Motion 


2.1. EQUATIONS OF LONGITUDINAL-SYMMETRIC MOTION REFERRED TO 
FIx—ED AXES 


Consider the motion of the aircraft with respect to the fixed rectangular axes 
defined in the Notation. Then, from Fig. 1, resolving with respect to x and z. 
we have 


mx=Pcos(6+y)—Lsin8-Dcosf, . ‘ (1) 

Gx i 


Moreover, with the conventional notation for derivatives, the elements of P, L, D, M 
for an infinitesimal increment dt of the time are 


dP=P.dz+P,dV, 

dL=L,dz+LydV + Lid6, 
dD=D,d2+DydV + Did6, 
dM = M.dz+M,dV +(Mo+K,) d6+K,d0. 


(5) 


But, as discussed in Appendix I, the aerodynamic terms in 4 may be neglected in 
the present analysis, stabilisation being then maintained purely through the control 
terms K,6 and K,6. Hence, the aerodynamic components of force and moment 
reduce to the instantaneous static terms while, for P, the characteristics of the 
turbo-jet engine indicate that, for the small ranges of z and V with which we shall 
be concerned, P, and Py may be regarded as constants. Hence, equations (1), (2) 
and (3) become 


mx =(P,+ P.z+PyV) cos (8+) —40p,V?S (C, sin B+ Cp cos), (6) 
mi =(P,+P.z+PyV) sin (6+ y)+40p,V2S . (7) 


where P, is the initial (sea level) value of P and V?= x? +Z?*. 


*Where gusts are involved, the value of V, under the conditions of symmetry considered here, 
must clearly be replaced by V?=(it+u)?+(z+w), u and w being, respectively, the horizontal 
and vertical velocity components of the gust. Equally, @ is then given by tan B=(z + w)/(x+u). 
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2.2. HOVERING STABILITY. BEHAVIOUR IN A VERTICAL GUST WITH No LAG 
OF THE CONTROL 


The simultaneous equations (6), (7), (8) can only be solved numerically in the 
general case, but when the aircraft is hovering the analysis is appreciably simplified 
for small displacements of the order which accord with satisfactory stability charac- 
teristics. The problem then reduces to one of ordinary mechanical or electrical 
circuit theory. At the same time, as already noted, and as will be seen in due 
course, the solution for this particular case reveals all the essential features pertinent 
to the more complex motion relating to non-stationary flight. Accordingly, the 
behaviour of the hovering aeroplane is first studied when subject to certain 
disturbances under the foregoing conditions. 


Consider, then, the effect of a gust. Examination of the relative influence of 
the horizontal and vertical velocity components shows that, for equal strength*, the 
vertical component has a predominant influence and, to shorten the analysis as much 
as possible, attention will be confined to the effects of a pure up gust or down gust. 
Further, to preserve the two-dimensional character of the problem, only the gust 
which, at a given instant, is everywhere of uniform intensity will be treated and, in 
accordance with generally accepted practice, its growth and decay will be 
represented by a sinusoidal function of the time. Relative to the gust velocity, it 
also appears (see Section 4.1) that the linear velocity components of the aircraft 
induced by the disturbance are sufficiently small to be neglected in so far as the 
aerodynamic forces and moments are concerned. Hence, the resultant air flow at 
any instant is assumed to be in the vertical direction while, for simplicity, the body 
axis is taken to be initially horizontal and the thrust axis vertical, in which case 


—a=B=y=2/2. 


Hence, for a small angular displacement @ with respect to the position of 
equilibrium, we have, since L then tends to zero and P to mg 


. . . . . @& 


9 being assumed sufficiently small to write sin? ~ 0, cos@~1. In accordance with 
the preceding observations, we also write 


w=w,sinot . ‘ . 12) 


in the interval O0<t<T,, where w, is the maximum gust velocity, 7, the duration 
of the gust and w the constant «/T,. 


*As is well known, atmospheric turbulence tends to be isotropic, i.e. the horizontal and vertical 
velocity components of a gust are approximately equal in magnitude, and the corresponding 
frequency distributions are essentially the same. 
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Then equation (11) may be written as 


= K, K, TP 


In terms of the Laplace transform 
co 

F(p)= 


equation (13) in the interval 0<t< T, takes the form 


20k, 


p(p? + 40”) 


(p+ 2ap+5*) f +(p+2a)0,+8,, 


where 9,, 9, are, respectively, the initial values of 8, 6, and, since we shall be 
concerned exclusively with the condition s?>a’, we write b?=s*-—a’. Then, 
factorising the quadratic on the left hand side of the last equation, 


f(p)= 207K, (p+2a)9,+9, 
p(p+a+ib)(p+a-—ib) (p+ 2iw)(p—2iw) (p+at+ib)(p+a-—ib)° 


(14) 


But, when t=0, 0=6=0. Hence, the second term of (14) is zero, and the inverse 
of (14) is then 


t 


I (a? — b? + — 2iab) 


+ conjugate } dt+ 


] 
+ (a? +b? — 4u? + 4iaw) + conjugate (15) 
0 


where “conjugate” indicates a similar term with the sign of i changed. 


Finally, writing = — b? + 40? =2a? — + 4w?, 


=a? + b?—4u*=s*— 


equation (15) yields 


6=e—" (C, cos bt + C, sin bt) + C, cos 2ot +C, sin 2ot+C,, . (16) 
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_ (A? +24) 

+ 4a"b*) (@ + 

v*K, 

2awK, 


Again, when t>T,, equation (12), and hence the first term of equation (14), are zero, 
while the second term of the same equation no longer vanishes, since, for continuity, 
6,, 8, now take the values of 6, @ at t=7,. Thus, 


= (p+2a)0,+9, 
(p+a+ib)(p+a-—ib)’ 


which yields 
sind (¢-7,)+0, 008 . (18) 


while, for continuity with (16) when t=T7,, 


6, =e-*"s[(aC, — bC..) cos bT, + (aC, + bC,) sin bT,] - 


6,=e-"e(C, cos bT, + C, sin bT,) + C, cos + C, sin 20T, +C,, 
(19) 
(C, sin 20T, —C, cos 2wT;). 


Hence, the solution of equation (11) for #, when w assumes the form of equation (12), 
is given by equations (16), (17), (18) and (19). The corresponding solutions of 
equations (9) and (10) follow immediately. Thus, in the case of equation (9), 
we have 


where C,, C, are the constants of integration. Then, for the interval O0<t< T, we 
obtain, on substituting (16) in (20) and (21), 


[ ac, +bC,) cos bt + (aC, — bC,) sin br | sin 2wt +£4cos 2wt—C,t+C,, 


a+b? 2 
(22) 
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2abC,} cos bt + {(a?— b?)C, —2abC,} sin br | + 
+ cos Qwt+ sin 2wt — CP + 


and, when t=0, x=x=0, in which case 


Cy 
a’ +b? 2w 
(24) 
(a’—b*)C,+2abC, C, 
(a? + b?)? 
Again, with respect to equation (10),z may be written as 
= = —k, Sin’ wf, 
in the interval O0<t<T,, where x,=op,w,?SCp/(2W). Hence, 
sin 
{ +D,t+D,, 
and when t=0, z=Z=0, therefore 
D,=0, D,= 
so that (1 —cos 2wf) — | (25) 


Clearly, equations (9) and (10) are not adequate approximations when t > T,, 
for then both x and z increase linearly with time, apart from a small oscillation of 
x which will be insignificant unless a—>0. In reality, however, the displacements 
will slowly decay due to the small aerodynamic resistance which has been neglected, 
and, equally, the corresponding damping in pitch will modify the solution (18) to 
equation (11). This, nevertheless, will only be significant if a is vanishingly small, 
when the solution given leads to an oscillation of constant amplitude instead of one 
which is slightly damped. After the gust has subsided, however, the motion is of 
no particular interest, for it is only the conditions prevalent in the interval 
0<t<T, which are characteristic of the problem under investigation, and it is 
then that the stability or instability of the aircraft will be decided. In the remainder 
of the discussion relating to pure gusts, therefore, attention will be confined to the 
interval 0<t< T;,. 
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2.3. HOVERING STABILITY. EFFECT OF CONTROL LAG IN A VERTICAL GUST 


The effect of lag in the control system could have been included in Section 2.2, 
but to avoid unnecessary elaboration it will be dealt with separately in the present 
section. Suppose that the time lags of the 6 and 6 responses may be expressed 
serially. Then, provided that the lag is a small proportion of the period of 
oscillation of the aeroplane, we may neglect all but first order terms and so replace 
6 by (6—«,6) and 6 by (6—<«,6), ¢, being the lag of the 6 system and «, that of the 
6 system. Hence, equation (13) becomes 


sin? wf, 


or 6+2a,6 +5,70=k, sin? of, . (26) 

_ 2a—s*e, 

where 24, = 
(27) 

Ky 
Thus, writing 


A?=a,? —b,?+ 4w?=2a,’ +40’, 
b,?—4w?=s,? — 


equation (16) has the modified form 
cos b,t+C,’ sin b,t)+C,’ cos 2ot+C,’ sin 2wt+C,’, . (28) 


for which equations (17) hold when C, is replaced by C,’, C, by C,’ etc., and a, is 
written for a, b, for b, x, for x,, A, for A and v, for v. The equations corresponding 
to (18), (19), (22), (23), (24) then follow immediately, as is obvious, by making these 
substitutions. It is also evident that, on the basis of the present treatment, lag of 
the system introduces no change in the form of the differential equation defining the 
longitudinal stability, i.e. equations (13) and (26) only differ in the magnitudes of 
the constant coefficients. It is, therefore, immaterial, for the purpose of general 
discussion, whether lag is considered to be present or not, and accordingly, attention 
is confined to equation (13) in the development which follows. 


2.4. HOVERING STABILITY. CHARACTERISTICS OF CONTROLS RESPONSIVE 
TO 8 or 9 EXCLUSIVELY 


It is of some interest to consider automatic controls which respond either to 
9 or to 6 only, since the resulting analysis shows clearly the essential contributions 
made by each of these elements to a control which will be satisfactory in a general 
manner. The former condition is merely a particular case (a=0) of the solutions 
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obtained in Sections 2.2, 2.3, so that no further consideration is necessary in this 
instance. On the other hand, the control purely responsive to @ requires separate 
treatment. Equation (13) then reduces to 


or, on integrating once, to 
t sin 2ot 
6+ 2a0=k, 40 (30) 


the constant being zero since 2=@=0 when t=0. Hence, solving this first order 
linear differential equation to satisfy the initial condition 6=0, t=0, we obtain, 
finally, over the interval O0<t< T,, 


=7 + cos 2ut Sin 2wt (31) 

Again, after the gust has subsided, i.e. for t>T,, equation (29) becomes 
6+208=0 . . 


or 6 + 2a0 =constant say, 
so that, on integration of (33), 


constant 


(34) 


Then © and the constant of integration in equation (34) may be determined from the 
conditions for continuity, namely, by substituting the relations for 6 and 6 given by 
(30) and (31) when t=T7, in (33) and (34). It is then found that 


where 2,, 2, are constants having the values 


_ sin 20T, 
(36) 


[1(7,- { sin 207, } 


Finally, in the interval 0<t<T7,, we have from equations (20), (21) and (31), 


= = -t) t+ a( 1 -cos 2ut- sin ) (1 ] 
(37) 
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08 (ato) } (38) 


in which the constants of integration are determined to satisfy the initial conditions 
x=X=0 when ¢t=0. 


2.5. HOVERING STABILITY. DISTURBANCES IN THE FORM OF A SUCCESSION 
OF SHARP PULSES 


Another type of disturbance of some practical interest which can be dealt with 
quite simply when the resultant air speed is zero, is one which persists with more or 
less regular frequency. For example, if the control jets are arranged in groups so 
disposed to give the requisite stabilising couples, it is possible that certain jets in a 
group may fail to operate correctly in that they introduce disturbing moments of a 
cyclical character through some fault in the servo-mechanism. Again, although the 
main jet-lift pipes may be arranged to give zero net thrust moment with respect to 
the centre of gravity of the aeroplane when the engines are functioning normally, 
any uneven running of one or more engines will tend to have a similar effect. 
Further, if moderate pulses of this kind occur with resonance frequency, they are 
liable to produce undesirably large displacements in a relatively short time should 
the automatic control be weak, but not necessarily inadequate, to meet other con- 
ditions. Inevitably such disturbances will have a very complex spectrum, but it is 
assumed here that the fundamental term is predominant, and can be represented by 
a periodic function which, in view of the essentially spasmodic nature of the 
disturbance, will be considered to have a discontinuous wave form. 


Under the present conditions, equations (6), (7) and (8) reduce to 


where P may now be a function of the time and AM is the disturbing moment 
about the axis, normal to the xz-plane, through the centre of gravity. In addition, 
as already seen in Section 2.2 on the evidence given in Section 4.1, the linear 
displacements x, z and their derivatives are of second order significance compared 
with the corresponding terms in pitch, so that essentially equation (41) remains. 
Then, in accordance with the foregoing remarks regarding the nature of the 
disturbance, we write 
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where H(t) is the unit step function, .=AM/B and 2g is the excitation period. 
Thus, writing for the moment, ¢ (¢) for the function on the right hand side of (42), 
we have 


(D? + 2aD +s’) 0=9 (0), 
or {D+(a+ib)} {D+(a—ib)} (0. (43) 


For convenience, let 6,(t) represent the complementary function and 6,(¢) the 
particular integral of equation (43). Then, 


6,()=e-“(Acosbt+Bsin bt), . . (44) 


with A, B as arbitrary constants, while, from the general relation for the particular [| 
integral of a second order inhomogeneous ordinary linear differential equation 


4 
co 
(a—ib)t r (a—ib)t_ a(a—ib)rq = 
{e 1} (t ra) {e }] 


—(a—tb)t 


(a—ib) 


pe — (a+b) t 


(a+ ib) (t) { i} rt 23 —1)'H (t-rq) { g(atid) ra } ] 


which, on reduction, yields finally 


6,()= EO (asin bt+ bcos on) } + 


+23 [asin b + bcos ra] | 
(45) 


Then 6=86, (t)+ 9, (0). . (46) 


But, when t=0, H(t)=H (t—rg)=0. Hence, the particular integral and its deriva- 
tives are all zero at t=0. Thus, @ is then equal to the complementary function 
6, (0), and, since 6=6=0 when t=0, it is seen from equation (44) and its first 
derivative that 


A=B=0. 


Therefore, 6=6,(1), . . (47) 


as given by equation (45). 
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3. The Case of Non-Stationary Flight 


3.1. NUMERICAL SOLUTIONS OF EQUATIONS (6), (7) AND (8) 


We will now consider the general condition when the aircraft has vertical and 
horizontal velocity components other than the values arising from a disturbance. In 
this case, the aerodynamic forces and couples become indeterminable functions of 
the time, and the simultaneous equations (6), (7) and (8) can only be solved 
approximately step-by-step. Take, first, equation (8). This may be written as 


. . . (48) 


where F(t) defines the quantity M/B. Then equation (48) has the solution 


[ A, 0s bt+B, sin bt+ F (7) e* sin b (t—7) dr | 


and by the rule for the differentiation of a definite integral 


t 


{a| F (7) e” sin b (t—7) dr— 


b 


6=-—e-* [ oA, +aB,) sin bt+(aA,—bB,)cos bt+ 


t 


-b F e* cos b (t—7) dr }] 


0 
A,, B, being the constants of integration. But, when t=0, 0=6=0. 


Therefore A,=B,=0, and hence, 


t 

—at 

F (7) sin b (t—7) dr, 


t 


t 
—at 
= [sin bt F (7) e* cos brdr —cos bt F (7) e” sin brdr | . (49) 
0 


the latter form for the solution of # being more convenient in the present discussion. 
Suppose now that over the time interval 0<t<t,, k equal sub-intervals r are 
taken, so that t,=kr. Also, let r be so small that, to a sufficient approximation, 
F(t) may be regarded as constant during the time lapse r and equal to its mean 
value. Thus, for (j-1)r<t<jr, write 


FO) =F tin] 


as the mean value of F(t) at the mean time 4(2j—1)r, where j=1, 2, 3, .. . k. 
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Then, substituting F(t) in (49), integrating over a sub-interval and summing over 
the total interval O0<(t<t,, we obtain, on writing t;_,=(j/—1)r, t; =jr, 


—at, 
sin bts, iF F@l, _{e*s(@cos bt, +b sin bt))— 


b(@ +b?) (a? + b?) 
— ej-1(acos bt,_, + b sin bt;_,)} — 


— Cos bt. = [F @) {es (asin bt;— b cosbt,) — e*i-1 (asin bt;_, — b cos bt;_,)} ] , (50) 
=1 -1 


from which it will be seen that, if t, is equated to the period 27/b, 


BA [F {ews (asin bt;—b cos bt;)—e*i-1 (asin bt;_,—b cos bt;_;)}, 
(51) 

while, when K,=a=0, there is the further simplification 


Equally, when t,=7/b, 6, is again given by equation (51), but the sign of the 
expression must be changed. The same condition applies in the case of equation 
(52). As will be seen in Section 3.2, where the numerical procedure is discussed, 
equations (51) and (52) can be used with great advantage to lighten the labour 
of computation. 


Now consider equations (6) and (7). In the interval 0<t < t,*, let x and z be 
represented by the quartics 


1 2 1 3 
Z=2,+2,t+ got? + Sat ‘ . (54) 


Then, by successive differentiation of (53) and (54) it follows that 


(55) 


80= Zo» 2:=(2) 8.=(Z), 


where the suffix 0 denotes an initial value, i.e. the value of a quantity at the beginning 
of a step, and given, therefore, by the end value of the preceding step. 


*It will be clear that, in a succession of such intervals, the origin is taken, for convenience, 
to shift from step to step. 
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Fig. 2. 
Experimental lift, drag and pitching moment data in a down gust. 


Z 


Further, tan B= sin B= cos B= (56) 
and + 2°, (57) 

while the aerodynamic coefficients 
(2), Cyo=f, (2), (2), (58) 


where 2=6—8 are assumed to be given empirically. Hence, through the first and 
second derivatives of equations (53) and (54), together with equations (56), (57) and 


(58), we obtain, on substitution in equations (6) and (7), two relations for h;, g, 


which may be expressed formally in terms of the given initial conditions as 


Unfortunately, equations (59) and (60) prove to be too involved for direct 
solution, but the values of h;, g, at each step may be obtained quite rapidly by 
successive approximation, provided that the aircraft is stable. This point is con- 
sidered further in the following section. 
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TIME t:SEC. 
Fig. 3. 


(A) a=0, b=2 sec.-', T,=3 sec. 
(B) a=0, b=2 sec.—', T,=5 sec. 
(C) a=1°7 sec.—!, b=0, T,=3 sec. 
(D) a=0'5 sec.-1, b=2 sec.—?, T,=3 sec. : 


3.2. NUMERICAL PROCEDURE. METHOD OF MEAN DEVIATIONS 


Given the initial values of x, z, @, together with the first derivatives of x, z, and 
taking P,, Py, f, (2), f. (@), fs (2) to be known empirically, the step-by-step solution of 
equations (6), (7) and (8) may be obtained as follows. Values of h, and g, in the 
interval concerned are first assumed*. Then x, z at the end of the interval follow 
from (53) and (54) respectively, and x, z by differentiation of the series, while the 
corresponding values of 8 and V are given by (56) and (57). Next assume a value 
of 2, and so determine C,, Cp, Cn. Now divide the step into a suitable number of 
equal sub-intervals r, calculate F () for each sub-interval and thence a first approxi- 
mation for @ from (50), or (51) and (52) when the conditions permit (see later). 
Hence, a second approximation for « follows from the relation 2=6—£. Finally, 
evaluate X, Z at the end of the interval by substitution in (6) and (7), and so, by 
double differentiation of (53) and (54), obtain second approximations for h,, g;. 
The calculations must then be repeated with the improved values of h,, g,, « until 
the approximations coincide to a satisfactory order of accuracy. 


*As the calculations develop, first approximations may be obtained by extrapolation from the 
data established, but in the early stages, when this is not a very reliable procedure, it is 
sufficient to start a new step with h,, g, assigned the values obtained finally in the 
preceding step. 
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When the control has a strongly damped element, as, for example, in Fig. 3 
(a=0°5 sec.~?, b=2-0 sec.-"), or in Fig. 12 (a=0-55 sec.-', b=1-80 sec.-'), the 
most convenient choice of interval is the natural period of oscillation 27/b, or some 
integral multiple of it, since the simplified equation (51) may then be used to 
calculate 6. If, however, there is an appreciable tendency to oscillation, as in the 
undamped, but stable, condition (e.g. see Fig. 3, a=0, b=2-0 sec.~-* and Fig. 12, 
a=0, b=1-80 sec.~'), the interval has to be reduced considerably, so greatly adding 
to the labour involved. The difficulty may be circumvented, however, by the 
following procedure* which results in a substantial economy of computation without 
loss of accuracy. 


Suppose that an approximate solution to the problem is first obtained by 
ignoring the oscillations in #, and that any convenient integral multiple of the natural 
period is again taken as interval. Then, the values of x, z and F(t) so obtained 
are strictly valid for a distribution of @ given, for example, by the lower envelope to 
the curve of Fig. 12 when a=0. Equally the upper envelope to @ is determined by 
calculating the half-periodic values. Both these estimates may be performed rapidly, 
since the much simplified relation (52) then holds. Finally, the mean deviation @,,, as 
indicated by the chain curve in Fig. 12, is used to re-calculate F (1), x, z as before, 
the process being continued until successive approximations agree with the requisite 
accuracy. In point of fact, it is found that the variations in 6, 8 and V, upon which 
the aerodynamic pitching moment, and hence F (1), depend, so combine as to make 
the second approximation to F(t) in very close agreement with the first. Hence, in 
practice, it is unlikely that the calculations need be repeated more than once. 


This method has also been carefully checked over one period against sample 
calculations which follow the actual oscillations of @ (the steps in this case were 
reduced to 1/10 of the period with sub-intervals of 1/60 of the period), and the 
comparisons, even for a severe gust, proved to be in agreement with an accuracy 
which is probably better than the general accuracy of the step-by-step method. 


This procedure of working on the mean deviation of #, when the motion is 
strongly oscillatory, involves little more time and labour than does the simpler 
calculation in the case of a control giving appreciable damping. 


4. Discussion of the Control Characteristics and the Conditions for 
Stabilisation 


4.1. GENERAL PROPERTIES OF 6 AND 6 CONTROLS 

Some examples illustrating the characteristics peculiar to the type of automatic 
control under consideration are now presented, and an examination is made of the 
conditions necessary for adequate stabilisation of the aeroplane when the normal 
aerodynamic damping becomes negligible. The results are divided into two groups. 


*Although the procedure has been found to give satisfactory convergence in the present 
calculations, and accords well with the checks referred to at the end of the section, no 
investigation has been made as to its general suitability and validity in step-by-step calculations 
of an oscillatory character. 
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The first set (see Figs. 3 to 8 inclusive; also Table I) refer to the hovering case, and 
the second (Figs. 9 to 14 inclusive) to some of the more interesting points arising 
from step-by-step calculations of the take-off flight path in a specific instance. 


As will be clear from Sections 2.2 to 2.5, the motion in hovering flight is 
proportional to the coefficients «,, «,, » which connect the mass and inertial charac- 
teristics of the aircraft with the forces and moments generated aerodynamically, or 
through some source of the kind discussed in Section 2.5. The calculations, there- 
fore, under such conditions, may be expressed quite generally, and accordingly, the 
relevant results, as presented in Figs. 3 to 8 and Table I, are independent of the 
geometry and mass distribution of the aircraft, and of the magnitude of the 
disturbance. On the other hand, when the aircraft is subject to a continuous 
displacement in addition to one of an oscillatory character, the aerodynamic loads, 
as already observed, then become dependent variables, and in these circumstances, 
the motion can only be analysed step-by-step for specific initial conditions. Those 
relating to Figs. 9 to 14 refer to an aircraft considered tentatively as a means of 
conducting some full-scale experiments on the type of stability problem discussed 
herein. The essential data were as follows: weight of the aircraft at take-off 
2,600 Ib.; corresponding wing loading 60 Ib./ft.?; sea level thrust (assumed to act, 
except where otherwise stated, at an angle y of 72:5°) 2,850 lb.; and thrust variations 
P., Py of, respectively, —0-0720 lb. ft.-', valid for 0<z< 5,000 ft., and 
— 1-325 Ib. sec. ft.-*, tenable over the range 0 < V < 150 ft./sec. 


Figure 2 also gives the lift, drag and quarter-chord pitching moment data upon 
which all the step-by-step calculations were based. These measurements do not refer 
to the aircraft in question, but were taken from some unpublished records of wind 
tunnel experiments on a model of similar form at a mean chord Reynolds number 
of 0:6 x 10°. The curves shown represent the average values of the coefficients from 
observations with (a) the normal wing, (b) the same wing with the span reduced by 
25 per cent. In each case the sets of coefficients differed very little when based on 
the respective wing areas. 


The first point to which attention will be drawn relates to the characteristics of 
the 6 and 6 components, and to the ways in which they are important for 
satisfactory stabilisation. Fig. 3 shows the oscillations in pitch under hovering 
conditions as the result of a vertical gust. Curves (A) and (B) are for a pure @ 
control, curve (C) for a 8 control only and curve (D) for a combination of the two. 
Ignoring, for the moment, the effect of lag, a control responsive to @ alone will, 
provided that it is of sufficient strength, preserve stability, but the aircraft continues 
to oscillate after the disturbance has ceased. Owing to the neglect of air resistance, 
this residual motion appears as a continuous oscillation of constant amplitude, 
whereas, in reality, it will decay, but, at the low velocities concerned, the damping 
will be so gradual that it will not materially alter the tendency of the aircraft to 
hunt. Fig. 12 also illustrates this characteristic very well during take-off in gusty 
conditions with the same form of control. Again, an adequate 6 control provides 
stability during a disturbance, but in contrast to the @ control, which preserves the 
original attitude of equilibrium, the 9 control causes the aircraft to assume a new 
Steady state, so that, unless correction is applied by manual intervention or by a 
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monitoring system, there will be a gradual divergence in, for example, a succession 
of gusts. The combination of 6 and @ response clearly avoids these disadvantages, 
and the less disturbed behaviour of an aircraft with the double form of control, as 
compared with that depending on @ only, during flight through a number of down 
gusts is well illustrated in Fig. 12. 


Of further interest is the way an aircraft reacts to gusts which persist for 
varying lengths of time. From Section 2.2, it is easily seen that, with zero damping 
(a=0), the effect of a gust will be a maximum when the duration of the gust is 
equal to the natural period of oscillation, i.e. when T,=27/b, or o=5/2. Thus, a 
dominant gust exists for every value of b. Examples demonstrating this feature 
are given in Fig. 4, from which it will be noted that, relative to the most severe 
condition, there is a much more pronounced variation of response to gusts of short 
existence than there is with long period gusts. Other calculations when a=1 sec.~? 
display essentially the same result, namely, a dominant gust which, for any value 
of b within the limits indicated in Fig. 4, has about the same duration as that 
corresponding to zero damping. Equally, there is, again, a greater variation of 
response to sharp gusts than to those of a more prolonged nature. It appears, 
therefore, that the aircraft exhibits a general uniformity of behaviour to gusts of 
various lengths over the ranges of a and b here considered, and that the value of T, 
for the dominant gust is virtually a function of b only. 


Figs. 5, 6, 7 present general curves of the displacements @n.x, x (T,) and z(T,), 
the first being the maximum value attained in the dominant gust* and the other two 
the values at the end of a gust. Estimates of «,, x, based on data from Ref. 1 for 
w, and Ref. 2 for B, and on the results of Fig. 2 for Cp, Cm, indicate that «, is not 
often likely to exceed about 0-1 sec~?, while «,, which is non-dimensional, also has 
a maximum limit in the region of 0-1. Thus, for instance, when b=1°5 sec.—', the 
displacements of the aircraft in the dominant gust, for which T,=4-2 sec., are: 
4°3°, x (T,) = 6°4 ft., z(T,)= 13-5 ft. when a=0, and 6,,.=1-7°, x(T,)=3°1 ft., 
z(T,)=13-5 ft. when a=1 sec.—', i.e. the average values of x, z during the gust 
period range from 1 to 3 ft./sec. Hence, since w, appropriate to the above example 
is of the order of 50 ft./sec., the assumption made in Section 2.2 that the components 
x and Z are aerodynamically insignificant, given adequate control, is quite justified. 


4.2. EFFECT OF LAG 


The crucial point from this aspect of the investigation is the absolute necessity 
for some degree of 6 control if stability is to be preserved when lag is present. This 
is evident from the first of equations (27) which shows that there is a critical value of 
«, for which the stability is neutral, i.e. for a,=0 or a=4s"<,, so that when the 


response to @ is absent (a=0), the aircraft is inevitably unstable however small 
the lag. 


*As already stated, the value of T, for the dominant gust proves to be virtually independent 

of a, and may, therefore, be taken as equal to the value 2x/b appropriate to a=0. In this 
case, however, equations (16) and (23) for 6, x respectively become indeterminate, since A, 
v are then both zero. Nevertheless, as Fig. 4 shows, the variation of 9 with T, for values 
close to 27/b is trivial, and the calculations were based on the gust times indicated by the 
curve AA of Fig. 4. 
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Thus, for any values of a and b, we have, since s?=a’ + b’, 


2. 
(crit)= (61) 
and further, for neutral stability, 
a’ +b? 
(62) 


Hence, given a and 5, the critical lag and the corresponding behaviour of the aircraft 
to a known disturbance may be determined according to the analysis of Section 2.3. 
Some examples based on the reasonable assumption that ¢,=¢,=¢, say, are recorded 
in Fig. 8. For a typical degree of lag, the minimum value of @ to give neutral 
stability is not large, and increases roughly in proportion to both « and b. Equally, 
the corresponding ratio of the displacement @,,.x with lag to that for zero lag is very 
nearly a linear function of ¢, but is not markedly dependent on b at the higher 
values. In the normal range of b, say 1-5 to 2:0 sec.—', and for a lag of 0°5 sec., 
Omax in the limiting case of neutral stability will be about 50-65 per cent. larger than 
the value when lag is absent. This is, perhaps, an extreme condition, but it 
illustrates the need for adequate damping from the control. 


In the further discussion, lag will again be ignored. 
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TABLE I 
(a) RESONANCE. LIMITING CONDITIONS FOR AN AMPLITUDE IN PITCH NOT EXCEEDING +10° 
a=0, b=1-047 a=0, b=2:094 a=0, b=3-142 
No. of Time No. of Time No. of Time 
half- lapse halj- lapse half- lapse 
cycles rq cycles rq sec,—? cycles rq sec.—? 
r sec. r sec. r sec, 
1 3 0:0955 2 3 0°1910 3 a 0°2865 
2 6 0:0478 4 6 0:0956 6 6 0°1434 
a 9 0:0318 6 9 0:0636 9 9 0:0954 
4 12 0:0239 8 12 0:0478 32 12 0:0717 
5 15 0:0191 10 15 0:0382 15 15 0:0573 
6 18 0:0159 12 18 0:0318 18 18 0:0477 
4 21 0:0136 14 21 0:0272 21 21 0-0408 
8 24 0:0119 16 24 0:0238 24 24 0:0357 
9 27 0:0106 18 at 0:0212 27 27 0:0318 
10 30 0:0096 20 30 0:0192 30 30 0:0288 
(b) DAMPED CONDITION. AMPLITUDE IN DEGREES AFTER A LAPSE OF 10 SEC. 
b=1-:047 b=2:094 b=3-142 
a=0°5 | a=1°0 a=0°5 | a=1°0 a=0°5 | a=1-0 
sec.-2 sec.—? 


0:0955 6°40 2°89 0°1910 6°57 3-20 0°2865 6°60 3°26 
0:0478 3-20 1°45 0:0956 3-29 1:60 0°1434 3-30 1:63 
0:0318 213 0:96 00636 2°19 1:07 00954 2°19 1-09 
0:0239 1:60 0°72 0:0478 1-65 0:80 0:0717 1:65 0°82 
0-0191 1:28 0°58 90382 1:32 0°64 0:0573 1-32 0°65 
0:0159 1:07 0°48 0:0318 1:10 0°53 0:0477 1:10 0°54 
00136 0°91 0°41 0:0272 0:94 0°46 0-0408 0:94 0°47 
00119 0°80 0:36 0:0238 0°82 0°40 0:0357 0°82 0°41 
0:0106 0°71 0°32 0-0212 0°73 0°37 0:0318 0°73 0°37 
0:0096 0°65 0°29 0:0192 0°66 0°32 0:0285 0°66 0-32 


Note: a and b are given in radian measure and have the dimension sec.-?. 


4.3. STABILITY DURING A SUCCESSION OF SHARP PULSES 

The influence of a repetitive disturbance, as discussed in Section 2-5, also 
demonstrates the importance of providing damping (a0) under certain conditions, 
namely, when the frequency of excitation approaches the natural frequency. 


In Table I(a), the allowable disturbing moment per unit moment of inertia for 
a pure @ control (system with no damping), such that the angle of pitch shall 
not exceed +10° after a specified number of oscillations, is recorded over a 
range of b. Inversely, Table I(b) gives the amplitude, with various degrees of 
damping, after a lapse of 10 sec. from the onset of the disturbance. Thus, whereas, 
for example, the undamped angle of pitch reaches 10° after 3 sec. when »=0°191 
sec.-? and b=2-094 sec.—'*, the corresponding damped displacement at the same 
values of » and b is only 6-6° after 10 sec. when a=0°5 sec.-’, or 3:2° when 
@=1-0 sec.-?. 


*The values of b given in Table I were adopted for convenience in order to make the excitation 
period 2g=2z/b a whole number. 
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4.4. STaBILIry FEATURES PECULIAR TO THE CASE WHEN THE DISPLACEMENTS 
ARE Not PURELY OSCILLATORY 


So far only the hovering condition has been discussed, when @, x, z are purely 
oscillatory, and x, Z may normally be ignored compared with the gust velocity w, 
but there are certain features of general interest when x, z exhibit a continuous 
positive or negative variation throughout the motion, e.g. during take-off or landing. 
As to the displacement 6, it has been noted that this aspect of the motion is much 
the same whether the aircraft is hovering or moving under its own power (cf. Figs. 3, 
12), and it is very useful that the simple hovering condition leads to results in this 
respect of wide application, since the stability of the aircraft (always, of course, with 
the restriction that the resultant air speed is sufficiently small to make the aero- 
dynamic damping negligible) is primarily dependent on the pitching motion 
generated. What happens with inadequate automatic control is portrayed in Figs. 9, 
10, 11. Thus, in Fig. 9 take-off flight paths have been plotted for the undamped 
condition (a=0) with, respectively, values of b equal to 0-90, 1:80 and 2-52 sec.-?*. 
In the first instance, the motion is clearly unstable. Due to the inadequate control 
moment, the pitch angle increases rapidly. As a consequence, the thrust inclination 
approaches the vertical (see Fig. 1), and then slopes in the opposite direction. 


*Again, as in the calculations relating to Table I, values of a and b for the step-by-step 
calculations were so chosen that the period, and therefore the interval, were whole numbers. 
This much simplifies the numerical work. Also, a was made zero, in order to shorten the 
calculations, whenever the immediate point under discussion was not significantly dependent 
on a. 
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Correspondingly, there is a decrease to zero of the horizontal thrust component, 
followed subsequently by a rapid increase in the reverse sense, while the vertical 
component remains very nearly constant at first, diminishes to zero as the body axis 
becomes normal to the flight path, and finally changes sign. It will be seen, therefore, 
that the initial forward and upward motion is steadily converted into one in a rear- 
ward and downward direction, so that the flight path takes the form of a loop, as 
indicated, with the aircraft in roughly the inverted position at the top of the trajectory. 
When, however, the control is adequate, @ varies very little (see, further, Fig. 10). 
Consequently, the horizontal and vertical thrust components remain of the same sign, 
and the flight path takes the desired form. Nevertheless, although the difference of 6 
between the controls with b=1-80 and b=2-52 sec.~* is very small, it is evident 
from Fig. 9 that the corresponding changes of flight path are much larger. This is 
due principally to the sensitivity of the horizontal thrust component to small 
variations of 6. It is also obvious that, in the unstable case (b=0-90 sec.~"), there 
will be a maximum for 2, as shown in Fig. 10, in the region where the aircraft body 
axis has become roughly normal to the direction of motion. Beyond this point the 
incidence will decrease again as the aircraft becomes inverted. At about the same 
time, the resultant air speed increases rapidly (see Fig. 11), since the thrust then acts 
principally along the flight path. Finally, had the calculations been extended to 
the culminating stage (an undertaking which would have involved a very con- 
siderable increase of labour owing to the rapid changes taking place), it is clear that 
they would have indicated a strong downward curvature of the flight path under the 
combination of gravity and inverse vertical thrust. 


The close accordance of the pitch characteristics, whether the aircraft is hover- 
ing or in general motion, has already been noted from Figs. 3 and 12. Equally, 
the earlier observation with respect to the small effect of a gust on x and z is again 
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apparent when translation of the aircraft is considered. To illustrate this point, an 
example is given in Fig. 13 of the flight path in calm and gusty air* when a=0, 
b=1'80 sec.-'. The close proximity of the curves agrees very well with the earlier 
conclusion derived from the analysis of the simple hovering case. Fig. 12 describes 
the corresponding pitching motion of the aircraft under the gust conditions. Of 
additional interest is the effect of adding some damping (a=0-55 sec.—") to the 
control. Here, the appreciable diminution of the oscillations in pitch (Fig. 12) 
results, for the reason given earlier, in a larger horizontal thrust component, and 
therefore the forward speed and displacement of the aircraft are increased relative 
to the corresponding undamped conditions. On the other hand, the motion 
vertically remains practically unchanged. Hence, the flight path becomes less 
steep. In general, therefore, any increase of damping or stiffness in the control 
leads to a more curved, i.e. a flatter, trajectory. 


Finally, in regard to this last feature, an interesting point is raised in Fig. 147, 
namely, the sensitivity of the initial flight direction 8 to the thrust angle y. This, 
of course, is again explained by the appreciable changes in the conditions hori- 
zontally for small differences of either # or y. Practically, it suggests that, where, 
for example, a take-off has to be made in a confined space, careful attention must 
be given to the setting of the jet pipes and ground angle of the aeroplane. Other 
discrepancies of possible significance, not here considered, may arise from lack of 
coincidence of the line of action of the jet thrust with the pipe axis (the most natural 
practical datum to assume) due to ground interference and asymmetry of the gas 
flow before discharge. 


*Triangular gusts were considered in contrast to the sinusoidal distribution assumed when 
hovering, since they are more convenient in the step-by-step calculations. Either form has 
substantially the same quantitative effect. 


tIn this case, it is quite adequate, for the purpose of demonstration, to assume infinite stiffness 
of the control (@ constant), when the calculations are particularly simple. 
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Appendix I 


ORDER OF THE TERMS IN @ IN EQuaTIONs (5) 


To show that the terms in # of equations (5) may be ignored under the conditions 
pertinent to the present investigation, it is noted first that either they are demonstrably 
negligible even under normal conditions of flight, or else that they arise predominantly 
from the action of the tailplane. To obtain the order of magnitude of the terms 
concerned, therefore, contributions from other parts of the aircraft will be neglected. 


Consider, now, the case when there is steady angular velocity 2 and 8 is constant. 
Then the tailplane experiences an increment of incidence 


Ind 


where a’ is the angle of incidence of the tailplane, J, is the distance of the tailplane 
lifting line from the centre of gravity and « is the downwash. 
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Further, by the usual assumption regarding downwash lag, namely that the 
instantaneous value of the downwash at the tailplane is the same as that appropriate to 
the lifting line of the wing at the instant when it occupied the spatial position now held 
by the tailplane, we have 


de 
| =} 
so that be’ = 1+ da (65) 


Again, when there is steady angular velocity B, with a constant, the increment of 
tailplane incidence is 


$a’ = > 


(66) 


so that when both 2 and £ are varying, we have, combining equations (65) and (66), 


where L’, D’, M’ are the lift, drag and pitching moment respectively contributed by the 
tailplane. If, further, we write 


L’, (Cy’, Cy’), 
M’= iop,V?Sc Cy’, 


where primes again denote measures related to the tailplane, neglect the downwash 
variation on account of the very high incidences involved and assume, for simplicity, 
that the tailplane chord is parallel to the wing chord, i.e. a’ becomes equal to a, we 
have for the instantaneous values of L, D, M 


(68) 


D ~ 40p,V°S(Co+ 


dC, 1,6 
M 4op,V28C (Cu+ 


The magnitudes of the various terms in equations (68) are now evaluated under the 
most severe conditions which can reasonably be expected. Whether they relate to a 
small but responsive aircraft, or to one of greater size and comparative sluggishness, the 
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value of [,6/V during a disturbance is unlikely to exceed 0-06. From tailplane experi- 
ments on the model referred to in Section 4.1 under the conditions there indicated, it 
also appears that over the interval 30°< a< 70° the coefficients C,’, Cy’, C,,’ are 
roughly linear functions of a with dC,’/da ~ +40-20 per rad.*, dCp’/da ~ 0-38 per rad. 
and dC,,’/da=—0-17 per rad. Accordingly, the following maxima are obtained for 


the 6 terms in equations (68) : 


aC, ib 
0023, 


Cy! 


dC 
+0-012, 


— 0-010. 


da V 


Relative to the maximum values of C,, Cp, C,, (Fig. 2), the corresponding terms above 
are all of the order of one per cent., so that even when the latter quantities are halved, 
which occurs when the incidence is about 80° for lift and 40° to 50° for drag and pitch, 
i.e. when considering virtually the full flight range of incidence, the error in neglecting 
the 6 terms is unlikely to be more than two per cent. 


Appendix II 
GyROSCOPIC YAWING COUPLE DUE TO PITCHING OF THE AIRCRAFT 


When the aircraft pitches, a yawing couple N, may be produced by revolving parts i 
such as the impeller, or impellers, of the turbo-jet power system. This effect depends [| 
upon the disposition of the various units. In general, if OX, OY, OZ are any set of [> 
rectangular axes with OX as the relevant axis of rotation, then, when the system as a 
whole experiences an angular velocity about OY, there will be a gyroscopic couple with : 
respect to OZ. Hence, in the case of the pitch and yaw interaction, the moment will be | 
principally generated by those units for which OX is substantially parallel to, but not 
necessarily in, the plane of symmetry, while, if N, is measured about the normal to the i 
longitudinal axis of the aircraft, only those revolving elements for which OX is parallel, } 
or approximately parallel, to the longitudinal axis will be significantt. Accordingly, in 
the following brief analysis, which is merely concerned with the order of the effect, 
secondary terms are ignored and the problem is regarded as essentially determined by 
the foregoing approximations relating to the disposition of reference axes and to the 
selection of significant mechanical units. Then, the magnitude of the couple N, with ‘ 
respect to OZ is 


where p is the angular velocity of the rotating parts (assumed constant and of the same 
magnitude and sense for each component), and / is the total polar moment of inertia. 
On the assumption, therefore, that the disturbance in pitch arises from a two-dimensional! 
vertical gust, as discussed in Section 2.2, and further, that any subsequent lateral motion 
may be treated independently of the longitudinal motion (an argument which is 
reasonably valid during the short interval of a disturbance), we have, when the 


*The lift curve has a sharp maximum at about %=50°, such that, for 30°< 4 < 50’, 
dC,’ /da~0-20 per rad., and for 50°< a<70°, dC,’/da—0°20 per rad. 


tNormally, the longitudinal axis will not be appreciably displaced from the horizontal, so 
that N, then refers roughly to the vertical, or OZ, axis of the fixed reference frame adopted 
in the main part of the paper. 
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Yaw due to gyroscopic action. 


yawing motion due to N, is reacted by automatic controls similar to those discussed in 
the case of pitch, the companion equation to (8), namely, 


where C is the moment of inertia of the aircraft about OZ, y is the angle of yaw, 
K,’, K,/ the control stiffness coefficients proportional to y and y respectively and 


x,’ 


N, _ 
26’ 


« 


Then the solution to (70) follows immediately from (49) simply by substituting y for 
6, a’ for a, b for b, F’ (t) for F (#) and the value of 6 in (69) by the differentiation of (16). 
The result may be expressed as 

t 
[ sin b’t e*’'Q(r) cos b’rdr — cos b’t | e*’Q(r) sin | 


cy 


pik, 


where 


2 
(r)=- {(aC, — bC,) cos br +(bC, + aC, sin br} — —=(C, sin 2wr — C, cos 2w7). 
Ky 1 


Fig. 15 is a plot of Ctmax/(plx,) as a function of b and b’ when a=a’=0 and 1< 6}, 
b’ <3, Wnex being the maximum angle of yaw in the gust. 
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value of 1,6 /V during a disturbance is unlikely to exceed 0-06. From tailplane experi- 
ments on the model referred to in Section 4.1 under the conditions there indicated, it 
also appears that over the interval 30°< a< 70° the coefficients C,’, C,’, C,,’ are 
roughly linear functions of a with dC,’/da ~+40-20 per rad.*, dCp’/da ~ 0-38 per rad. 
and dC,,’/da=~—0-17 per rad. Accordingly, the following maxima are obtained for 


the 6 terms in equations (68): 


dC,,’ 

0-010. 
Relative to the maximum values of C,, Cy, C,, (Fig. 2), the corresponding terms above 
are all of the order of one per cent., so that even when the latter quantities are halved, 
which occurs when the incidence is about 80° for lift and 40° to 50° for drag and pitch, 
i.e. when considering virtually the full flight range of incidence, the error in neglecting 
the 6 terms is unlikely to be more than two per cent. 
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GyYROSCOPIC YAWING COUPLE DUE TO PITCHING OF THE AIRCRAFT 


When the aircraft pitches, a yawing couple N, may be produced by revolving parts 
such as the impeller, or impellers, of the turbo-jet power system. This effect depends 
upon the disposition of the various units. In general, if OX, OY, OZ are any set of 
rectangular axes with OX as the relevant axis of rotation, then, when the system as a 
whole experiences an angular velocity about OY, there will be a gyroscopic couple with 
respect to OZ. Hence, in the case of the pitch and yaw interaction, the moment will be 
principally generated by those units for which OX is substantially parallel to, but not 
necessarily in, the plane of symmetry, while, if N, is measured about the normal to the 
longitudinal axis of the aircraft, only those revolving elements for which OX is parallel, 
or approximately parallel, to the longitudinal axis will be significantt. Accordingly, in 
the following brief analysis, which is merely concerned with the order of the effect, 
secondary terms are ignored and the problem is regarded as essentially determined by 
the foregoing approximations relating to the disposition of reference axes and to the 
selection of significant mechanical units. Then, the magnitude of the couple N, with 
respect to OZ is 


where p is the angular velocity of the rotating parts (assumed constant and of the same 
magnitude and sense for each component), and / is the total polar moment of inertia. 
On the assumption, therefore, that the disturbance in pitch arises from a two-dimensional 
vertical gust, as discussed in Section 2.2, and further, that any subsequent lateral motion 
may be treated independently of the longitudinal motion (an argument which is 
reasonably valid during the short interval of a disturbance), we have, when the 


*The lift curve has a sharp maximum at about %=50°, such that, for 30°< 4 < 50°, 
dC,’ /da=0-20 per rad., and for 50°< a<:70°, dC,’/dae—0°20 per rad. 


tNormally, the longitudinal axis will not be appreciably displaced from the horizontal, so 
that N, then refers roughly to the vertical, or OZ, axis of the fixed reference frame adopted 
in the main part of the paper. 
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“< 
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O5 20 25 3-0 
COEFFICIENT b:SEC' 
Fig. 15. 


Yaw due to gyroscopic action. 


yawing motion due to N, is reacted by automatic controls similar to those discussed in 
the case of pitch, the companion equation to (8), namely, 


where C is the moment of inertia of the aircraft about OZ, y is the angle of yaw, 
K,’, K,/ the control stiffness coefficients proportional to y and y respectively and 


K,’ 


36° 


Then the solution to (70) follows immediately from (49) simply by substituting ¥ for 
6, a’ for a, b for b, F’ (t) for F (#) and the value of @ in (69) by the differentiation of (16). 
The result may be expressed as 


t 
—a’t 

[ sin b’t cos b’rdr — cos b’t | e*Q(r) sin b’rdr | 

6 

where 
{(aC, — bC,) cos br + (bC, + aC, sin br} — (C, sin 2w7r C, cos 2w7). 
Ky 1 


Fig. 15 is a plot of Cym.x/(plx,) as a function of b and b’ when a=a’=0 and 1 <5, 
b’ <3, Wax being the maximum angle of yaw in the gust. 
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Ww. S. COLEMAN 


From an examination of the data on current turbo-jet engines, it would appear 
that the quantity pl/P varies between about 1 and 2 rad.ft.sec. Further, taking the 
weight W as roughly equal to the thrust P, and the radius of gyration in yaw from the 
data of Ref. 2, it is estimated that pI /C will not normally exceed 1-5 rad. sec.~*. Hence, 
with the maximum value assigned to x,, namely 0-1 sec.~? (see Section 4.1), a figure of 
0-15 rad. sec-* is obtained for the upper limit to plx,/C. Thus, from Fig. 15, it will 
be seen that at roughly the lower limits for b and b’, say, b=b’=1°5 sec.—', the value 
of Wmex is not likely to exceed about 5°, and when a £0 it will be even smaller. Hence, 
the assumption of two-dimensional motion advanced in Section 1 appears to be 
reasonably justified. 


Errata—May 1955 (Vol. VI, Part 2) 


In the paper by E. W. GraHaM, B. J. BEANE and R. M. LicHer, “ The Drag of Non- 
Planar Thickness Distributions in Supersonic Flow,” the authors have asked for the 
following addition to be made on p. 111:— 


“In Figure 7, the equivalent Sears-Haack body of revolution should be of fine- 
ness ratio 5. The central body of revolution should be approximately double-conical 
in shape (i.e. two cones placed base to base), and the maximum thickness of the 
central body, for the given wing geometry, is approximately half the maximum 
thickness of the Sears-Haack body.” 


In the paper by G. N. Lance, “ The Lift of Twisted and Cambered Wings,” the 
following corrections are necessary : — 


p. 154: The argument of the cosh in (15) should read 3 | ; 


p: 155; Insert a minus sign on the right hand side of (18). 


p. 156: The left hand sides of the equations on lines 9, 10 and 11 and 13 should 
read b,/a,,, b,/a,,, C./%,,, and +7 respectively. 
The right hand side of (25) should be — 1, 


pris: Replace (27) and (28) with 


ab,’ [2a?K —(1 +a?) E] + ab,’ [a?K (Sa?— 3) + 2E(3 Sa +.a*)] + 
+a’c,/[2a*K . (27) 


ab,’ [aK (5 — 3a”) + 2E (3a* — Sa? + 1)] + ab,’ [2a*K — (a? + a*) EE] + 
+a°*c,’ [a*K (3 a?)—2E]=0 (28) 


The left hand sides of (30) and (31) should read —xb,’ and +b,’ 
respectively. 


p. 161: In the expression for ¢,,, the first «,, should be preceded by a minus sign 
and the term —2rn? should read +217. 


The author wishes to express his thanks to Mr. Julian Kainer and Mr. Hitch for 
drawing his attention to these errors. 


The Aeronautical Quarterly 


328 


ear 
the 
the 
ice, 

of 

ill 
lue 
ice, 

be a 
on- 
the 
ne- 
cal 
the 
um 
the 
uld 
27) 
28) 

2 
ign 
for 
erly 


i 
| 
| : : 
| 4 
| 
| 
| > j 
| 
| 
| 
} 
| 
| 
| 
| 
‘ 
| ) 
| 
| 
} 
| 
} q 
| a 
| 
| 
| 
| rs 
> 
= 
| 
| 4 
| 
| 
4 
} 
| 
: 
| 


INDEX TO 


AERODYNAMICS 

See also Channel Flow, Compressor 
Blades, Photography, Stability, Sting 
Supports. 

Note on a Null Method for the Direct 
Measurement of Pressure Coefficients, 
L. H. Tanner, p. 241 (November 1955). 

On an Approximate Method of Calculat- 
ing Pressures on Non-lifting Head 
Shapes at Supersonic Speeds, H. K. 
Zienkiewicz, p. 31 (February 1955). 

On the Unsteady Motion of a Slender 
Body through a Compressible Fluid, 
L. E. Fraenkel, p. 59 (February 1955). 
Errata, p. 164 (May 1955). 

Supersonic Flow Past Slender Bodies with 
Discontinuous Profile Slope, L. E. 
Fraenkel and H. Portnoy, p. 114 (May 
1955). 

The Distribution of Pressure on an Aero- 
foil, in a Stream with a Spanwise 
Velocity Gradient, W. A. Mair, p. 1 
(February 1955). 

The Drag of Non-Planar Thickness Distri- 
butions in Supersonic Flow, E. W. 
Graham, B. J. Beane and R. M. Licher, 
p. 99 (May 1955). 

Addendum, p. 240 (August 1955). 
Erratum, p. 328 (November 1955). 

The Lift of Twisted and Cambered Wings 
in Supersonic Flow, G. N. Lance, p. 149 
(May 1955). 

Erratum, p. 328 (November 1955). 


BEANE, B. J., E. W. GRAHAM and 
R. M, LICHER 


The Drag of Non-Planar Thickness Distri- 
= in Supersonic Flow, p. 99 (May 


Addendum, p. 240 (August 1955). 
Erratum, p. 328 (November 1955). 


BISHOP, R. E. D. and D. C. JOHNSON 
Some Properties of Nodes in Vibrating 
Systems, p. 46 (February 1955). 
CHANNEL FLOW 
Compressible Subsonic Flow in Two- 
dimensional Channels, L.C. Woods. 
Part I: Basic Mathematical Theory, 
p. 205 (August 1955), 
Part II: The ae of the Theory 
to Problems of Channel Flow, p. 254 
(November 1955). 


COLEMAN, W. S. 


Stabilisation of the Aeroplane in Flight at 
Zero or Very Small Air Speeds, p. 295 
(November 1955), 

COMPRESSOR BLADES 

Some Effects of Mechanical Damping on 
the Self-Induced Aeroelastic Oscillations 
of an Axial Compressor Blade, I. M. 
Davidson, p. 81 (May 1955). 


VOLUME VI 


DAVIDSON, I. M. 
Some Effects of Mechanical Damping on 
the Self-Induced Aeroelastic Oscillations 
of an Axial Compressor Blade, p. 81 
(May 1955). 


DONALDSON, I. S. 
The Effect of Sting Supports on the Base 
Pressure of a Blunt-Based Body in a 
Supersonic Stream, p. 221 (Aug. 1955). 


ELECTRICAL ANALOGUE 
See Structures. 


ELLIOTT, D. E. 
The Measurement of High Speed Air 
Velocity and Temperature using Sound 


Wave Photography, p. 181 (August 
1955). 


FATIGUE 
See also Photoelasticity. 


Effect of Yielding on the Fatigue Proper- 
ties of Test Pieces containing Stress 
Concentrations, K. Gunn, p. 277 
(November 1955). 


FRAENKEL, L. E. 


On the Unsteady Motion of a Slender 
Body through a Compressible Fluid, 
p. 59 (February 1955). 

Errata, p. 164 (May 1955). 


FRAENKEL, L. E. and H. PORTNOY 
Supersonic Flow Past Slender Bodies with 
Discontinuous Profile Slope, p. 114 
(May 1955). 


GODFREY, D. E. R. 
Normal Loading on a Wedge-Shaped 
Plate, p. 196 (August 1955). 


GRAHAM, E. W., B. J. BEANE and 
R. M. LICHER 


The Drag of Non-Planar Thickness 
Distributions in Supersonic Flow, p. 99 
(May 1955). 

Addendum, p. 240 (August 1955). 
Erratum, p. 328 (November 1955). 


GUNN, K. 
Effect of Yielding on the Fatigue Proper- 
ties of Test Pieces containing Stress 
_— p. 277 (November 


HOLISTER, G. S., H. T. JESSOP and 
C. SNELL 


Photoelastic Investigation in Connection 
with the Fatigue Strength of Bolted 
Joints, p. 230 (August 1955). 


INTERPLANETARY FLIGHT 


Dynamic Problems of Interplanetary 
Flight, Derek F. Lawden, p. 165 
(August 1955). 


M4 

| 

329 


JESSOP, H. T., C. SNELL and 
G. S. HOLISTER 
Photoelastic Investigation in Connection 
with the Fatigue Strength of Bolted 
Joints, p. 230 (August 1955). 


JOHNSON, D. C. and R. E. D. BISHOP 


Some Properties of Nodes in Vibrating 
Systems, p. 46 (February 1955). 


LANCE, G. N, 
The Lift of Twisted and Cambered Wings 
in Supersonic Flow, p. 149 (May 1955). 
Erratum, p. 328 (November 1955). 


LAWDEN, DEREK F, 


Dynamic Problems of Interplanetary 
Flight, p. 165 (August 1955). 


LICHER, R. M., E. W. GRAHAM and 
B. J. BEANE 


The Drag of Non-Planar ‘Thickness 
Distributions in Supersonic Flow, p. 99 
(May 1955). 

Addendum, p. 240 (August 1955). 
Erratum, p. 328 (November 1955). 


MAIR, W. A. 
The Distribution of Pressure on an Aero- 
foil, in a Stream with a Spanwise 
Velocity Gradient, p. 1 (February 1955). 


PALMER, P. J. and S. C. REDSHAW 
Experiments with an Electrical Analogue 
or the Extension and Flexure of Flat 
Plates, p. 13 (February 1955). 


PHOTOELASTICITY 
Photoelastic Investigation in Connection 
with the Fatigue Strength of Bolted 
Joints, H. T. Jessop, C. Snell and G. S. 
Holister, p. 230 (August 1955). 


PHOTOGRAPHY 
The Measurement of High Speed Air 
Velocity and Temperature using Sound 
Wave Photography, D. E. Elliott, p. 181 
(August 1955). 


PORTNOY, H. and L. E. FRAENKEL 


Supersonic Flow Past Slender Bodies with 
Discontinuous Profile Slope, p. 114 
(May 1955). 


REDSHAW, S. C. and P. J. PALMER 


Experiments with an Electrical Analogue 
for the Extension and Flexure of Flat 
Plates, p. 13 (February 1955). 


REVIEWS 


Linearized Theory of Steady High-Speed 
Flow, by G. N. Ward. Reviewed by 
W. A. Mair, p. 240 (August 1955). 

Studies in Mathematics and Mechanics, 

resented to Richard von Mises, edited 
y Garrett Birkhoff, Gustav Kuerti and 
Gabor Szegé. Reviewed by M. J. Light- 
hill, p. 164 (May 1955). 


ROBINSON, J. R. 
The Buckling and Bending of Orthotropic 
Sandwich Panels With All Edges 
Simply-Supported, p. 125 (May 1955). 


SANDWICH PANELS 
The Buckling and Bending of Orthotropic 
Sandwich Panels With All Edges 
Simply-Supported, J. R. Robinson, 
p. 125 (May 1955). 


SNELL, C., H. T. JESSOP and 
G. 8S. HOLISTER 


Photoelastic Investigation in Connection 
with the Fatigue Strength of Bolted 
Joints, p. 230 (August 1955). 


STABILITY 


Stabilisation of the Aeroplane in Flight at 
Zero or Very Small Air Speeds, W. S. 
Coleman, p. 295 (November 1955). 


STING SUPPORTS 
The Effect of Sting Supports on the Base 
Pressure of a Blunt-Based Body in a 
Supersonic Stream, I. S. Donaldson, 
p. 221 (August 1955). 


STRUCTURES 

See also Fatigue, Photoelasticity, Sand- 
wich Panels. 

Experiments with an Electrical Analogue 
for the Extension and Flexure of Flat 
Plates, P. J. Palmer and S. C. Redshaw, 
p. 13 (February 1955). 

Normal Loading on a Wedge-Shaped 
Plate, D. E. R. Godfrey, p. 196 (August 
1955). 


TANNER, L. H. 
Note on a Null Method for the Direct 
Measurement of Pressure Coefficients, 
p. 241 (November 1955). 


VIBRATION 
See also Compressor Blades 
Some Properties of Nodes in Vibrating 
Systems, R. E. D. Bishop and D. C. 
Johnson, p. 46 (February 1955). 


WOODS, L. C. 


Compressible Subsonic Flow in Two- 
dimensional Channels. 
Part I: Basic Mathematical Theory, 
p. 205 (August 1955). 
Part II: The ——— of the Theory 
to Problems of Channel Flow, p. 254 
(November 1955). 


ZIENKIEWICZ, H. K. 


On an Approximate Method of Calculat- 
ing Pressures on Non-lifting Head 
Shapes at Supersonic Speeds, p. 31 
(February 1955). 


INDEX 
| 
330 


iges 
opic 
dges 
son, 
tion 
ited 
t at 
Base 
in a 
ison, 
and- 
ogue 
Flat 
haw, 
aped 
gust 
irect 
ents, 
ating 
eory, 
eory 
. 254 
ulat- 
Head 


